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COMPUTER PROGRAM FOR CALCULATION OF A GAS TEMPERATURE 
PROFILE BY INFRARED EMISSION - ABSORPTION SPECTROSCOPY 

by Donald R. Buchele 
Lewis Research Center 

SUMMARY 

A computer program to calculate the temperature profile of a flame 
or hot gas is presented in detai l. EmphasisJLs on profiles found in jet 
engine or rocket engine exhaust streams containing H 2 O or CO 2 radi- 
ating gases. The temperature profile is assumed axisymmetric with an 
assumed functional form controlled by two variable parameters. The 
parameters are calculated using measurements of gas radiation at two 
wavelengths in the infrared. The program also gives some information 
on the pressure profile. A method of selection of wavelengths is given 
that is likely to lead to an accurate determination of the parameters. 

The program is written in FORTRAN IV language and runs in less than 
60 seconds on a Univac 1100 computer, i 

INTRODUCTION 

Infrared emission and absorption measurements at two or more 
wavelengths provide a means of determining a gas temperature profile 
along a path through the gas by using a racjiation source and receiver 
located only at the ends of the path (ref. 1).< rln; the ease of a jet engine 
nr, rocket engine exhaust stream, the instrumentation may be outside of 
the stream. The method has been termed the ’‘spectral- scanning 
method. ” Practical instruments that permit use of this technique are 
described in references 1 through 4. , 

Tests of the method with heated pure gases in a furnace of 60 cm 
path length (ref. 2), and then with flames over a 22 cm path length 



2 


(ref. 3) gave good results. An analysis in reference 5 showed the effect 
of radiometric errors on the accuracy of the computed profile. The 
analysis assumed the partial pressure of the radiating gas was a con- 
stant over the entire path, with the constant determined from radicjmetric 
measurements. More recently, the feasibility of using only emission 
measurements and an estimated gas composition to determine industrial 
furnace temperature profiles over a 10-meter path length at atmospheric 
pressure has been shown (ref. 6). 

The computer program presented herein involves an iterative, 
cyclic computation in which an initial assumed (arbitrary but realistic) 
temperatui 0 profile is altered in shape until the computed emission and 
absorption that it would yield^agree as closely.as possible with the actual 
instrumental measurements of emission and absorption. The procedure 
requires the use of a large-scale computer. These procedures have been 
outlined in principle in previous publications; but this report presents an 
actual program that has been found to produce useful results. 

The present work also extends the work of reference 5 by providing 
greater latitude in the temperature- profile shape that can be treated and 
by relaxing restrictions on the shape of the associated pressure profile. 
Thus, whereas previous studies have principally treated profiles with a 
substantial temperature difference over the optical path, the present 
study also includes the equally important case of an almost isothermal 
profile; and whereas the work of reference 5 has assumed constant 
pressure over the path, the present work allows an additional degree of 
freedom in the description of the pressure profile. 

The program presented provides a means of determining temperature 
from experimental measurements of emission and absorption at two or 
more wavelengths. But it also provides a means, before any experiments 
are performed, of selecting those wavelengths whose use will lead to the 
best attainable accuracy. This prior selection is very important, because 
a poor selection of wavelengths may lead to unacceptably high uncertainty 
in the temperature and may even lead to an erroneous solution of the 
equations. The existence and avoidance of this incorrect solution for the 


temperature profile is demonstrated. Thus, this report prtmdes the 
potential user with a means of determining whether this method will be 
sufficiently accurate for his particular application, as well as providing 
a means of finding the solution if the L..ethod is adopted Symbols are 
listed in appendix A. 


MATHEMATICAL ANALYSIS 

Method of Determining Temperature and Pressure Profiles 

The method is iterative with each iterative step consisting of three 
principal parts: 

1. With current profiles of temperature and pressure, compute the 
gas radiance and transmittance at each wavelength of measurement, 

2. Record the difference between computed and measured gas radi- 
ance and transmittance at each Vvavelength of measurement. 

3. Correct the current profiles in (1) to reduce the differences in 
(2) to zero. 

Functional forms of the temperature profile T(x) and the partial 
pressure profile p(x) of the absorbing gas are assumed in advance. The 
temperature profile is expressed in terms of two unknown parameters,/ 
and the pressure profile in terms of one unknown parameter, A mint^ 
mum of one set of measurements, at each of two wavelengths, is required 
to yield a solution for the unknown parameters, each set actually in- 
volves four radiometric measurements (i.e. , four instrument readings). 
The computational program then determines, from the simultaneous 
solution of equations containing actual radiometric measurements, the 
values of these parameters that malce computed values of radiance and 
transmittance equal to measured values. 

Additional wavelengths of measurement may improve accuracy, in 
such case, a least-squares method of solution is employed. The results 
obtained are the values of the parameters when the ^'residual'*, which 
is the sum of the squares of the differences between measured and com- 
puted values of radiance, is a minimum. 
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The functions used in this program may be modified as desired for 
a particular application. However, tests during development of the 
pri*)gram have shown that it is impractical to characterize the tempera- 
ture profile by more than two unknown parameters because the accuracy 
of the solution becomes unacceptably low. 

The profiles . - In the present program, the profiles are taken to be 
axisymmetric, with shapes chosen to resemble the transverse pressure 
and temperature distributions often found in ducts. The functions are 
expressed in terms of values of the dependent variable at the axial, or 
centerline, position (subscript c) and at a wall position (subscript w), 
and the value of a shape parameter or for temperature or pres- 
sure, respectively. 

1. When the temperature profile is far from isothermal, it is taken 
as 


T(x) = T^-(Tg-T^)|y|^/’>T 


( 1 ) 


where 


y = (2x/L) - 1 0 X < L 

It is assumed that T,„ is measured independently. Parameters T„ 

w c 

and rirji are determined by the computational program. They are then 
termed the dependent parameters. This family of curves is shown in 
figure 1(a). 

2. When the profile is nearly isothermal, except for a tliin boundary 
layer at each end, the temperature profile is taken as 

T(x) = T^ - (T^ - T*)|y| ^ (2a) 


y = (2x/L) -1 XjSx^sL-Xj 


where 
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3 

and 

T(x) = + (Tj - T^)z + (mxj + - Tj^)z(z - 1) 0 :s z :s 1 (2b) 

where 


z = x/xj for 0 :< X Xj^ 

z = (L - x)/Xj^ for L - Xj :S X :s L 

Tj =T(xj)=T(L-Xj) 

m = (dT/dx)^^^^ = 2(Tg - Tj)/j^^(L - 2Xj^ 

This definition of T(x) provides continuity of both temperature and tem- 
perature gradient at x - x^ and at x = L - Xj^. This family of curves 
is shown in figure 1(b). Using this definition, 77 and Xj^ are preassigned 

some reasonable values and the computational program determines T* 

w 

and T- as the dependent parameters. It is assumed that T „ is meas- 

L W 

ured independently. 

3. The partial pressure profile of each absorbing gas is assumed 
to be of the form 

I I 

p(x) = P(. - (Pg - P^)ly| P (3, 


where 


y = (2x/L) - 1 0 :< X L 

The profile is shown in figure 1(c). Different initial values may be assumed 

for parameters p„, p„ . and ? 7 „ for each gas. The choice of th initial 
o >v p 

values of the parameters may be aided by a computation of the partial 
pressure of the reaction products. These pressures can be computed for 
a range of reaction product temperatures and a range of partial pressures 
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and reactant composition using the computer program of reference 7. 
From these data a choice for initial values of the pressure parameters 
can be made. (The program of ref. 7 may be obtained from its authors, 
on written request, by sending a magnetic tape at least 1200 ft long and 
by specifying the type of computer. Parts of the program not used in 
this application may be removed after the program has been trans- 
ferred from tape to cards. ) 

Tests during the development of the program for computing temper- 
ature and v'essure profiles have shown that only one dependent pres- 
sure parameter can be determined with acceptable accuracy. This is 
best done by designating the parameter p^ j of the principal absorbing 
gas (subscript 1) as an independent parameter. The wall pressure 
p,„ 1 is selected to be directly proportional to the center pressure, and 
the shape parameter n ^ is chosen judiciously and held constant. 
Assumed parameters for all other absorbing gases (subscript k) and 
all nonabsorbing gases (subscript j) determine the initial pressure pro- 
files [Pu(x)]j, and [PjCx)]^. 

Next, these profiles are made dependent on the profile P]^(x). The 
absorbing gases initially have profiles that determine ratios 
[Pj^(x)/p^(x)Jq. Then, as p^(x) is changed by the program, these other 
profiles are also clianged by the linear relation 

Pj(x)(pj^(x)/p^(x)]q (4) 

After the absorbing-gas pressure profiles are determined, the non- 
absorbing-gas pressure profiles Pj(2^) adjusted to make the sum of 
all gas pressures equal to the measured static pressure Pg which is 
assumed to be independent of x, so that 

Ps =DPkW-^X)PjW 

k j 

An additional constraint is to maintain the relative amounts of each non- 
absorbing gas as given by the initial profiles, that is, in the ratio 
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[Pj(.x)/^ Pj(x)3o 
j 

Thus the pressures of the nonabsorbing gases are given by 

= [Ps -D Pk^’^o [Pj(=‘)/D Pj(='Bo (5) 

k j 


The residuals . - Initial values of all parameters of the profile func- 
tions are tabulated, and the resultant sets of values of gas radiance 
and transmittiavje are computed for two or more ivavelengllis from 
radiometric equations. The compuiod quantities are denoted by 

Q JLj 

and The corresponding measured quantities at the same wave- 
lengths are denoted by and The difference between the 

computed and measured quantities is called the ^^residual. ” The re- 
sidual for gas radiance may be expressed as a fractional value 


Rn = ("Nl - (6) 

It is preferable to use gas absorptance 1 " '^L Pl^ce of trans- 

mittance Tj^ for the transmittance residual because a-^ is more nearly 
proportional to the gas pressure and to the gas absorption coefficient. 

The residual for gas absorptance is defined as 






(7) 


The solution . - The solution is reached by reducing the residuals to 
an acceptably small value by appropriately changing the two variable 
parameters in the equation for the temperature profile and the one vari- 
able parameter in the equation for the pressure profile. This is done 
iteratively using the Newton-Raphson method (ref. 8, p. 175) because 
the equations are nonlinear. For two wavelengths of measurement each 
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ilorativc step involves the solution of the following three simultaneous 
equations for the increments Ay. to the current value of the three vari- 
able parameters represented by y., 







?y4 


Ay. - -R 


Ai,N 



/ 



= -%N / 


( 8 ) 





Oi 


Note that only one wavelength is needed for R. to determine one 

A ^ j CJC 

pressure parameter. The derivatives in equation (8) are obtained numer- 
ically. The final value of the residuals is zero. 

There are cases where there are more measurements than there are 
variable parameters. For these cases, where the final values of the 
residuals are not necessarily zero, a damped least-squares modification 
of equation (8) is used to find a solution at which 



a 



- a minimum 


( 9 ) 


where there are as many residuals R^^ as there are i wavelengths, 
and there is only one residual R^. The damped least squares method 
is used in this program. Its derivation is in the second part of appendix B. 
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In addition to the solution for the profile parameters, an f>; 3 timate 
is made of the error in the parameters that would ’'e caused by stated 
random errors of radiometric measurements and tabulated absorption 
coefficients. (Typical orders of magnitude for these errors are 1 per- 
cent in radiation and 5 percent in absorption coefficient. ) The error 
equation is derived in appendix B. 


Measurement of Gas Radiance and Transmittance 

The actual gas radiance Nj^ and transmittance are derived 
from measurements of detector (radiometer) output voltage V that is 
assumed to be proportional to incoming radiant power. The factor of 
proportionality (calibration factor) of the detector includes the effect of 
all optics between the radiometer element and the gas, but it need not 
be known explicity if it remains constant for all measurements. The 
optics include, in addition to windows, lenses, and stops, a monocliro- 
mator of bandwidth A\ centered on wavelength X. The following 
measurements of detector output voltage V a’e made: 

(a) Vq, where a calibrated source of radi ice is placed at the 
gas location, in place of the gas. 

(b) Vj^, due to radiation from the gas. 

(c) V 2 , where a stable comparison source radiates through the 
gas path, but the gas is absent. 

(d) Vg, when the comparison source radiates through the gas path 
in the presence of the gas. 

By chopping the radiation from the comparison source at a unique, 
distinctive frequency, voltages Vj and V 3 are distinguished from each 
other, and voltage Vg ooes not include gas radiation (ref. 1 , p. 28). 

These four measurements yield 





V. 

V 


Tl = 


0 



1 


( 10 ) 
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Computation of Gas Radiance and Transmittance 

The gas radiance, integrated over a narrow wavelength band AX 
(tenter ed at X and over a path length L, is 



where is the Planck radiation function 


Ci Co/XT 

= -(e ^ 

X 



( 11 ) 


and T is the gas transmittance over the path from 0 to x. The de- 
tailed justification for equation (11) is presented in reference 5. 

The gas transmittance for a band AX is described by a band- model 
formulation for a nonisotherraal gas based on equation (5) of reference 9, 
Tins formulation has been justified and applied in reference 2, and sum- 
marized in reference 10. For profiles with a large variation in tem- 
perature, the formulation as modified in references 3 and 11 is adopted. 
Other band- models can be found in the literature (ref, 12). 

The transmittance r for a single gas is a function of the dimen- 
sionless optical depth u and is given by 

r 1 


T = exp 



(12) 


For a mixture of q absorbing gases the combined transmittance is the 
product 


T - 



1=1 


T, 


(13) 



The optical depth of a path length x is 


u 



( 14 ) 


v/here k is the absorption coefficient at local temperature T and par- 
tial pressure p of the absorbing gas and 17 is an interpolation param- 
eter. The parameter /3 is the local value of the line overlap parameter 
and is given by jS - 2ny/d where y is the average line width, and d is 
the average line spacing. The definitions of u and r given here differ 
slightly in form from those of references 3 and 11 for reasons of compu- 
tational expediency. 

A modification in references 3 and 11 replaces the constant average 
i3g of reference 2 of ref. 4) by a cumulative average produced by 
weighting /S(x) with the absorption coefficient, so that 



(15) 


The interpolation parameter 17 is given by one of two formulas 


u + Aj 

V = u > Ap. 

u + Ag 

> 

uAo 

7] = U Ac 

A5 + UA4 


(16) 


where 
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Values for A^, A 2 , Ag, A^, and Ag are listed in table I in terms ol 

The computation for the initially assumed profiles T{x) and p(x) 
first determines all values of k, (3, and /3 , as functions of x from 
(he band model data and equation (15). For example, tabulations of k 
and (^/d)^ for HgO, CO 2 , and CO may be obtained from reference 13 
or, for HgO, from references 14 and 15. (The subscript s here de- 
notes the value at standard temperature and pressure, ) Values of k, 
in the tables are denormalized to the local gas temperature and partial 
pressure by 

k=kp^ (18) 


where p is in atmospheres and T isinkelvin. Similarly, the tabu- 
lated values of (y/d)^ = ^/{2n) may be used to obtain the local value as 
it is affected by the line -broadening effect of all gas components in 
accordance with the formula 


/3 






(19) 


where p^, pj^ are the partial pressures, and 0 ;? , are the broaden- 
ing coefficients for the k absorbing gases and j nonabsorbing gases. 

For the absorbing gas, values of a* are given in reference 16 and are 
listed in table II. 

The path length 0 x ;< L is divided into increments Ax^ = x^ - x._. 
(i = 2, . . . ) and the values of k, /3//3 , /3 are determined at each 
station. Then, starting with the values r/^ •■= 0, u^ =0, and = 1 at 
X - 0, the following sequence of computations is made for each successive 
increment of path: 

(1) Find u^=u._^+AUj with equation (14), using k^, (/3//i^g)j. 
and 
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(2) Find 17^^ with equations (10) and (11), using u^. Linear inter- 
polation in table I is used to find to A^, since the computation is 
relatively insensitive to these quantities. Equation (11) is used to find 

^ 5 * 

(3) Find with equation (12) using u^ and (/3^)^. 

For the numerical computation, equations (14) and (15) are written 




1 



(14a) 


i 



(15a) 


Accuracy of the Solution 

There are several sources of error that affect the accuracy of the 
solution. First there is the basic limitation of the functional forms 
assumed for the temperature and pressure profiles. Obviously not all 
situations encountered in practice can be accommodated by these func- 
tions with acceptable accuracy. For those cases other functional forms 
must be used. This basic limitation is demonstrated in the section 
"Application of the Program to an Experiment" where an example is 
worked out for a temperature profile that cannot be matched exactly by 
the assumed functional form for the temperature. 

Another source of error is generated by the nature of the set of 
iteration equations for the parameters, equation (8). On the left side. 



the matrix of derivatives may be poorly conditioned. The conditioning 
of the matrix is affected by the selection of the wavelengths of meas- 
urement, Proper selection of the wavelengths is treated in the next 
section. 

The residuals on the right side of equation (8) are not exact quan- 

c c m m 

titles due to errors in their components Tj^, and Tj^. 

The measured quantities and contain an experimental 

random measurement error and the computed quantities and 

are in error because of uncertainties in the tabulated gas parameters 
/3_ and k_. There are also errors due to deficiencies in the band 
model used in the gas transmittance calculation. These deficiencies 
in the model are covered in the references cited in the previous section 
and will not be covered further here. Thus, even if the computer pro- 
gram succeeds in finding a solution that reduces the residual values to 
zero, the computed parameters are still in error due to the residual 
errors. By using equations (6) and (7) it can easily be shown that the 
residual errors at the solution, and 6R^ are given by 


6R 


N “ 







(20a) 


(20b) 


where and are the errors in 

c in 

Tj^, and Tj^, respectively. 

To proceed further it is necessary to specify the errors in the 
measured and computed quantities. It will be assumed that all meas- 
urements have the same probable error given by | 6^V/*^v| and that 
the computed values for Nj^ and have random errors due to the 
uncertainty |6k_/k„| in published k„ values. 

' o o ' S 
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The uncertainty in the value of /3g will be neglected. Thus, only 

two random errors will be considered, |6*^V/^V| and |Sk_/k„| . 

I ' s s 

Using equations (10) and the definition ^ 

easily shown that 



(21a) 


(21b) 


(21c) 


To estimate the random errors in the computed quantities, assume a 
simplified transmittance equation given by = exp(-kpL). It can 
then be shown that 



(21d) 


The above quantities when substituted into equations (20a) and (20b) 
give for the squares of the random residual errors 


6R 


2 

N " 




(22a) 


ORIGINAL 

OP POOR quality 
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The random residual error 6Rj^ (eq. (22a)) depends on the error 
6k/k. An alternative form that eliminates this dependence is to substi- 
tute for equation (6) the expression 


R 


N 




(23) 


where is an apparent radiance defined by 
By repeating the previous procedure it can be shown that 



2 2 

The quantity oRj^ is less than 6R^ whenever 


(24) 


(25) 


(in > 2(e^vrvf (26; 

(6k/k)2 

— o o 

Because experience has shown that 6Rjj. is less than 6R^ in most 
practical situations, the computer program presented herein uses Rj^ 
given by equation (23) instead of Rj^ given by equation (6) for the re- 
siduals in equations (8) or (9). 

Equations (25) and (22b) give a measure of the uncertainties in the 
zero values for the residuals at the solution. These expressions are 
used in the computer program following the procedure outlined in ap- 
pendix B to compute the random errors in the temperature and pressure 
profile parameters. 
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Selection of Wavelengths of Measurement 

The wavelengths differ principally in their associated gas absorp- 
tion coefficients. Figure 2 shows the absorption coefficient k for the 
2.7 jum band of K^O. At wavenumbers near the center of the band 
(3704 cm“^), k decreases with temperature; at the wings, k increases 
with temperature. Substantial differences in k exist, especially at a 
low temperature. Another difference due to wavelength is caused by a 
greater increase of the Planck- function radiance with temperature at 
the shorter wavelengths. Since the variation of k with wavelength is 
the phenomenon that enables the profile to be determined, large differ- 
ences would be expected to be favorable. 

The wavelengths are selected to maximize accuracy of the solu- 
tion, and equally important, to avoid an incorrect second solution that 
can exist. This preliminary selection is performed through an auxil- 
iary computer program WAVES. The selection procedure will be illus- 
trated with both modes of the temperature profile function, equa- 
tions (1) and (2). 

Avoiding an incorrect solution . - The WAVES program uses the 
following procedure: For assumed values of the two variable temper- 
ature parameters in equations (1) or (2), as applicable, and one variable 
pressure parameter in equation (3), the quantities and are 

computed for each wavelength. These computed quantities are substi- 
tuted for presumed error-free measurements and of the 

assumed profiles. Then, taking one wavelength at a time, one param- 
eter w of the temperature profile is assigned a value different from 
that assumed. The quantities and change. Then, new 

values of the other temperature parameter and the pressure parameter 

are found by iteration that restore the equalities = ^Nt and 

c in ^ 

= Tj^. This procedure is repeated for a number of assigned 

values of the temperature profile parameter w, and the results are 
tabulated by the program. 
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For each wavelength, the tabulation can be plotted as a curve with 
the independently varied parameter w as abscissa and the dependent 
temperature parameter as ordinate. The intersection of any two curves 
is a two- wavelength solution for a temperature profile. 

As an example, assume a path length of 100 cm, hydrogen -air 
combustion at 1 atm total pressure having H 2 O pressure parameters 
p = 0.35 atm, p„, = 0.30 atm, and n„ = 0. 5. Computations are made 
for wavelengths in the 2. 7 iim band of H 2 O. Results are plotted in 
figure 3 for the profile of figure 1(a) with n,p as the independent param- 
eter and T as dependent, and in figure 4 for the profile of figure 1(b) 
with as the independent parameter and as dependent, with no 
wall boundary layer (xj^ = 0). To facilitate plotting and to emphasize 
differences in the curves, the ordinate Y of the graph is tabulated by 
the program as a function 


Y = Z + ?(w-Wq) (27) 

where Z is the dependent parameter, w is the independent parameter, 

Wq is the originally- assumed value, and ^ is chosen by the program 

as described in appendix C. All curves intersect at the true solution. 

Some pairs of curves intersect at a second point that is an incorrect 

solution. For example, in figure 4, wavenumbers 3500 and 3100 cm"^ 

intersect at approximately T^ = 2200 K and T^ = 1840 K. A solution by 

the computer program converged to T* = 2223 K and T., = 1873 K from 

w c 

a starting point T^ = 2300 K and T^ = 1800 K. To assure convergence 
to the true solution, a pair of curves should have increasing separation 
with distance from the true solution. Accuracy is greatest when there 
is a large difference of slope between a pair of curves. Where two 
curves become parallel, accuracy approaches zero. 

Inspection of figures 3 and 4 show that no single pair of wavelengths 
is satisfactory for all profiles. Acceptable pairs of wavenumbers in 
figure 3 are (3100, 3550 cm”^) and (3100, 3775 cm"^). The 2000 K 
almost- isothermal profile that led to figure 4 shows an acceptable pair 
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(3500, 3100 cm"^) for T* below 2000 K, and an acceptable pair 

H W 

(4100, 3400 cm"-") for T* above 2000 K, 

w 

Tests of the pressure profile when the two pressure parameters 
p^ and p^ were assumed to be variables, and only was vari- 
able, gave curves nearly alike for all wavelengths. Thus, poor ac- 
curacy of the computed pressure profile makes this technique unac- 
ceptable; it is therefore necessary to limit oneself to only one variable 
pressure parameter. The reasons for this poor accuracy are (1) that 
radiance and absorption coefficient usually vary more strongly with 
temperature than they do with pressure and (2) that the measurement 
of Tj^, equation (10) which is the principal source of information about 
the pressure, is sensitive only to the integral of pressure over the 
path and not to local values of pressure, whereas Nj^ (which is the 
principal source of information about temperature) ^ sensitive to local 
temperature. 

Maximizing accuracy of the solution . - The program that led to 
figures 3 and 4 does not provide an estimate of the error for pairs of 
wavelengths. This estimate is provided by the main program PROFIL 
(MODE=l). The error estimate is made by pairing one principal 
wavelength with each of all the other wavelengths to be tested. For 
each dependent profile parameter, the program tabulates the error in 
the parameter due to assigned errors of the radiation- measurement 


voltages V and the tabulated absorption coefficient k . 

O 

As an example, such tabulated values are plotted in figure 5 for the 
parameter T^, with the profiles and some of the wavelengths of figure 3. 
The assumed errors were 1 percent in V and 5 percent in k^. Another 
example with the same assumed error is plotted in figure 6 for T 
with the profiles and some of the wavelengths of figure 4. 'These errors 
are for profiles with T = 2000 K, and the range of values covered by 
the abscissa parameter. The errors are consistent with the slope - 


difference of curve pairs in figures 3 and 4. The error is minimized 
in figure 6 with wavenumber pairs (4100, 3400 cm~^) above 1950 K, 


and (3500, 3100 cm" ) below 1850 K. 


OKIGlNAtPA^ 
OF POOR QUAUl.it 
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All four of these wavenumbers have been used simultaneously in a 
least squares solution by the program PROFIL (MODELS). The result- 
ing error also plotted in figure 6(b), shows improvement at all values 

of T* . The effect of path length on the error at T* = 1200 and 2000 K 
w w 

is listed in table III, for this four -wave length combination The error 
varies approximately as the inverse square root of the path length, in 
this example. 


PROGRAM DESCRIPTION 

The accuracy of the temperature profile found with measured data 
is dependent on the pressure and temperature profiles, path length, 
and wavelengths selected for measurement. It is desirable as a pre- 
liminary step to assume realistic values for the path length and pres- 
sure and temperature profiles for the absorbing gas, and then to select 
wavelengths and estimate the accuracy of the temperature profile 
parameters as was illustrated in the Selection of Wavelengths of Meas- 
urement section. 

A separate main program WAVES (described in ^’Avoiding an incor- 
rect solution”) is designed for the selection of wavelengths. Using these 
wavelengths, a principal main program PROFIL gives results in three 
modes of operation: 

MODE=l: This program computes the error in the dependent pro- 
file parameters that would be caused by assigned random errors of 
radiometric measurements and tabulated absorption coefficients. For 
given values of all profile parameters, gas radiance and transmittance 
are calculated and the parameter errors are computed as in appendix B. 
Actual radiometric measurements are not required. The estimate may 
be made for any pair of wavelengths consisting of one ’'principal” wave- 
length and any one of up to seven others. 

MODE=2: This program computes values for the dependent profile 
parameters using the iterative procedure described previously. Limits 
are assigned to the dependent parameters by the user; if a limit is 
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reachod durin^^ the iteration, that parameter becomes fixed, and the 
iteration proceeds for the remaining dependent parameterSo When a 
least- squares solution is used, the magnitude of the final root- mean- 
square residual R (eq. (9)) is also provided. In addition, using the 
.same procedure as in MODEsl, the error in the dependent parameters 
tliat would be caused by assigned random errors in the radiometric 
measurements and tabulated absorption coeflicients is computed 

MODELS: This program is used to find the change of temperatur(> 
profile that would be caused by a change in one of the independent pro- 
file parameters. A secona set of data cards is required that incor- 
porate the cliange, Starting with an assumed profile and gas-property 
data, the resulting radiance and transmittanc e that would be measured 
are calculated. The second set of data cards is then introduced to 
change one of the independent profile parameters, and the computation 
in MODE -^^2 is performed using the previous radiance and transmittance 
values to yield a revised profile. This process is useful In identifying 
independent profile parameters to which the temperature profile is 
unusually sensitive. 

The flowchart for program PROFIL is figure 7 After reading input 
data and calculating some initial conditions, the mam program calls 
subroutines TEMP that calculates profiles, and TRAMS that calculates 
gas transmittance and radiated flux. The subroutines are principally 
a computation using equations presented in the text. The equation 
numbers appear on comment cards in the listings. 

The test result DONE-1 is obtained when (1) the number of iterations 

reaches a preset number, or (2) when the root-mean-square sum of all 

residuals is less than 10~^, or (3) when the fractional change of the 

variable parameters per iterative step has become less than 10“^ for 

-3 

T^ and T^, and less than 10 for p^ and Urj,. Criterion (3 j is 
needed when there are more residuals than variable parameters, be- 
cause the residuals then reach a nonzero minimum, 

A subsidiary loop of moderate complexity determines the variable 
parameter increments and tests the parameters for exceeding their 
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limits. The loop contains a control index number KEY that manipulates 
the index of FORTRAN DO statements as follows. The variable profile 
parameters in equation (1) are n,p, T^,, p^, T^. They are represented 
by a FORTRAN subscripted variable PA(K) that is controlled by a DO 
index K=;K1, K2, K3. An index initially K=l, 3, 1 selects parameters 
ii;j„ T^, Pg for the profiles of figure 1(a) and K=2, 4, 1 selects param- 
eters T^, Pp, T^ for the profiles of figure 1(b). 

Any parameter in either profile can be fixed when it reac lies an 
assigned limit by changing the index parameters Kl, K2, K3 This is 
done by setting Kl, K2, K3 equal to the subscripted quantities 
KEYl(KEY), KEY2(KEY), KEY3(KEY) that can be manipulated. The 
subscript KEY can select any combination of numbers for Kl, K2, K3 
as shown in table IV. The variable parameters for each KEY are also 
listed. By designating the profiles of figure 1(a) as PROF=l and the 
profiles of figure 1(b) as PROF=2 the table is valid for PROF=l and 
PROF =2 when subsequent additions are made, Kl=Kl-l+PROF, and 
K2=K2-l+PROF. The initial profile with no limits reached is KEY=1. 

Before entering the loop, KEY=1, and Kl, K2, K3, are set as in 
the table. This frees the currently fixed parameters to return inside 
their limits during every iterative step. Thus, on first entering the 
loop, all three parameters (corresponding to KEY=1) are incremented 
and compared with their limits. H the K'th parameter is within its 
limit, a limit-indicating variable LIM (K4)=0 where K4=K-PROF+l 
(to be valid for PROF=l and PROF=2). If the KTh parameter is outside 
the limit, LIM(K4)=1. The value of KEY is a number from 1 to 8 
given by 


KEY=l+4* LIM(l)+2* LIM(2)+LIM(3) 

If this new value of KEY is not 1 or 8 the preceding procedure of this 
paragraph is repeated only once from the start of the loop with the 
new value of KEY, and the corresponding Kl, K2, K3 to find new 
parameter increments and tentative new parameter values. The loop 
is then exited, the tentative new parameter values are then accepted, 
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and tlic itoration proceeds. At the solution the error of the profile 
parameters is estimated. The prog’ram listing with a FORTRAN symbol 
list is appendix D. The subroutine lisnngs are appendixes E and F. 

The flow chart for program WAVES is figure 8, After reading the 
input data, the program (following the procedure described in Selection 
of Wavelengths of Measurement) first sets and ^ 

for the initial profile. Then fell': v> three «iested loops, An outer DO 
loop for selecting each given waveiengtln a DO loop for assigning given 
values to one independent variable profile* parameter: and an itcratii»n 
loop like PROFIL that finds the other two dependent variable* parameters. 
If a variable parameter reaches a limit during iteration, the* iteration 
is ended. After completion of the two DO loops the results are tabu- 
lated. The program listing, with FORTRAN symbols not used in pro- 
gram PROFIL, is appendix G. 

Successive problems in the same run may omit input data iliat is 
not changed, (Both PROFIL and WAVES programs, ) The READ data 
is omitted by not putting IREAD^l in the NAMELIST data , In the 
NAMELIST, the data-inpiit parameters of temperature, T^„, T_, and 

and pressure, p^, p^, and n^ are given starting values that are 
not changed by the program because they are replaced by different 
internal symbols for c'^mputation. After computation the values of the 
internal symbols are printed, but with the corresponding symbols of 
the initial data- input parameters, 

PROFIL PROGRAM USE 
Input 

The input consists of four types of data: (1) computer program c on- 
trol constants (2) measurements of gas radiance and transxnittance at 
selected wavelengths (3) tabulations from the literature of gas trans- 
mission parameters at selected wavelengths and (4) independently deter- 
mined or estimated profiles of partial pressures oi the gas reacUon 
products. 
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The collected input data are read in with a FORTRAN NAMELIST 
for single -subscript quantities, followed by a FORTRAN READ 
STATEMENT for multiple- subscript quantities. An example in 
appendix H shows this input data in the same format as punched on 
cards, for MODE=l, five wavenumbers, and two gases, H 2 O and N 2 . 

The data conform to the following specifications. 

NAMELIST DATA . - The following list gives the FORTRAN symbol 
with its subscripted maximum dimension and eiiUivalent algebraic 
symbol, in the NAMELIST sequence. 

FORTRAN symbo l Algebraic Description 

svrnbol 


MODE MODE=l finds error estimate 

of pairs of wavenumbers con- 
sisting of the first wavenumber 
paired with each succeeding 
wavenumber 

MODE -2 finds solution to given 
measurements, with error 
estimate 

MODE =3 finds gas radiance and 
transmittance for given input, 
then requires NAMELIST re- 
peated that may have changes 
from the first NAMELIST, and 
finds solution the same as 
MODE =2 

PROF PROF=l profile figure 1(a), var- 

iable parameters TC, PC, 

NT. PROF =2 profile fig- 
ure 1(b), variable parameters 
TW, TC, PC, with given bound- 
ary layer wall temperature TWB, 
thickness X(JB) 


FORTRAN symbol 

Algebraic 

symbol 

Description 

IREAD 


IREAD=1 requires READ data 
for GAS, ALPH, KTAB, BTAB. 
With IREAD=1 omitted, READ 
data must be omitted 

PATH 

L 

Test section path length, cm 

11 (8 maximum) 


Number of wavenumbers for com- 
putation (when IREAD=1, 11 
must equal number of data 
cards for KTAB and BTAB) 

WAVE (8) 

X 

_i 

Wavenumber cm 

MRL(8) 


Measured radiance W, cm" ^ 
cm~^ ster”^ (eq. (10)) 

MTAUL(8) 


Measured transmittance (eq. (10)) 

ERRV 


Percent error of radiation meas- 
urement 

ERRK 


Percent error of tabulated ab- 
sorption coefficient 

TW 

rii rp^ 

w’ w 

Gas temperature, K, at v^all. 
For PROF=2 with boundary 
layer, T^ approximates free 
stream temperature near the 
wall 

TC 

Tc 

Gas temperature, K, on axis 

NT 

Urp 

Gas temperature profile equation 
exponent 

TWB 

T 

w 

Gas temperature, K, at wall for 
PROF =2 with boundary layer 
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FORTRAN symbol 

Alg'^braic 

symbol 

Description 

JB (minimum is 1 at wall) 


Index for X(JB) giving distance 
X from wall where boundary 
layer ends, for PROF --2 

UPLIM(4) 


Upper limit of profile param- 
eters TW, TC, PC, NT 

BOTLIM(4) 


Lower limit of profile param- 
eters TW, TC, PC, NT 

LI (3 maximum) 


Number of absorbing gases 

L2 (6 maximum) 


Number of absorbing plus non- 
absorbing gases 

PS 

Ps 

Gas static pressure, constant 
along optical path 

PC(6) 

Pc 

Partial pressure of each gas, on 
axis 

NP(6) 


Pressure profile equation ex- 
ponent of each gas 

PW(6) 

Pw 

Partial pressure of each gas, 
at wall 

NK (8 maximum) 


Number of tabulated values of 
KTAB read in for data 

TKTAB(8) 


Temperatures at which values 
of KTAB are tabulated in 
READ data 

NB (8) maximum) 


Number of tabulated values of 
BTAB read in for data 

TBTAB(8) 


Temperatures at which values of 
BTAB are tabulated in READ 
data 
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FORTRAN symbol Algebraic Description 

symbol 

J1 (41 maximum) Number of stations along optical 

path. Must be odd. Omit if 
Jl=21 in DATA -statement of 
program is acceptable. 

Z(J) Station distance ratio, from wall 

to center, as fractional value 
0 to 1 . Omit if list in DATA- 
statement of program is 
acceptable. Maximum 
J=(Jl+l)/2 

ITER Maximum number of iterations 


READ statement da ta. - The following list gives the FORTRAN sym- 
bol with its subscripted maximum dimension, the format of each card, 
and the equivalent algebraic symbol. 


FORTRAN 

symbol 

GAS(6) 

ALPH(3, 6) 


Format Algebraic Description 

symbol 

6A6 Names of all gases 

(12F6. 0) Q!* Broadening factor for up to 3 ab- 

sorbing gases by up to 6 total 
number of gases. Card se- 
quence is one card for each LI 
absorbing gases. Each card 
lists the broadening factor for 
each of L2 gases. 


/ 


t 
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Description 


Absorption coefficient at S, T- P. 
for up to 11=8 wavenumbers, 

LI =3 gases, NK=8 tabulated 
points. Card sequence is one 
set of cards for each of LI 
absorbing gases. Each set has 
II cards corresponding to the 
wavenumbers in NAMELIST. 
Each card lists absorption coef- 
ficients corresponding to tem- 
peratures TKTAB in NAME- 
LIST. 

Line- broadening coefficient at 
S. T. P. for up to II =8 wavenum- 
bers, LI =3 gases, NB-=8 tabu- 
lated points. Card sequence 
same as for KTAB above. Each 
card lists broadening coeffi- 
cients corresponding to temper - 
ature TBTAB in NAMELIST. 

The two gas transmittance parameters KTAB and BTAB depend on 
temperature. For H 2 O and COg, they are tabulated in references 10 
to 12 with KTAB and BTAB normalized to S. T. P, The program inter- 
polates between tabulated values to find a tabulated value at local tem- 
perature, then denormalizes to obtain KK and BB at local temperature 
and pressure in accordance with equations (18) and (19). Linear inter- 
polation is used. If the profile temperature is outside the temperature 
range of the tabulation, KK and BB are limited to the extreme tabulated 
value. It is t-hus desirable to prepare the tabulation with a tcfmpe nature 
range running from the lower to upper limits appearing on olher input 
cards for the profile temperatures TC, TW, and TWB. 


FORTRAN 

symbol 


Format 


Algebraic 

symbol 


KTAB(8, 3,8) (12F6.0) k 


BTAB(8, 3,8) (12F6.0) 
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Output 


An example of the output shown in appendix H for MODE =1 and 
PROF =2 consists of four parts: 

1. The input NAMELIST and READ data. 

2. Estimated error of the profile parameters caused by the re- 
siduals, equations (‘jO) and (25), due to the errors of measurement 
ERRV and of absorption coefficient ERRK. The example has ERRV=1 
percent and ERRK=5 percent. The example shows four wavenumbers 
3100, 3400, 4100, and 3000 cm"^ successively paired with the first 

wavenumber 3500 cm”^. The first pair have residuals R^-r =2-0 and 

- 1 ^ -1 
3.2 percent at 3500 and 3100 cm , and R^ =4.4 percent at 3500 cm 

These residuals cause the listed three errors for each parameter, and 

their rms total. 

3. Transmittance at each wavenumber, 

4. Temperature and pressure profiles. 

An example shown in appendix I for MODE =3 was in the same run 
and followed the preceding example. The READ statement data and most 
items in the NAMELIST were not changed. Thus, only two data cards 
were required. The first four wavenumbers were used in a least squares 
solution. Appendix I shows the input data in the same format as punched 
on cards. The first NAMELIST has MODE=3 and 11=4. The required 
second NAMELIST has a change in the nonvariable pressure parameter, 
p^ = 0. 36 replacing p^ = 0. 30. The output consists of four parts. 

1. The first input NAMELIST and the second NAMELIST that differs 
by the new p^ = 0. 36 and an internal change to MODE =2. Because 
IREAD=1 was omitted, the READ data were not required or printed out, 

2. Each iteration step lists the parameters, the residuals Rj^ for 

each wavenumber, the residual R^ for the first wavenumber, therms 

total residual, and a step size that is the total fractional change of the 

variable parameters. In this example the 20 percent change of p„, 

from 0. 30 to 0. 36 caused errors 6T = 1. 2 K or 0. 06 percent, 6p = 

-0. 019 atm or 5,4 percent, and 6T = 14 K or 0.5 percent. Because 

w 
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four wavenumbers were used (11=4) the iteration was terminated by the 
small step size, as the residuals became stationary and remained 
nonzero. 

3. Estimated error of the profile parameters caused by the residuals 

and The percent errors are = 3, 5 percent, 6p^ = 

7. 1 percent and 5T „ = 6. 6 percent. These errors are much larger than 

w 

those listed in (2) above caused by the change of p„ » 

w 

4. Temperature and pressure profiles. 

The output in MODE=2 is the same as MODE =3 without the first NAME- 
LIST. 


WAVES PROGRAM USE 
Input 

The input includes all quantities listed for PROFIL. Three additional 
quantities in the NAMELIST select the independent temperature param- 
eter and assign values to it that differ from the initial value. The inde- 
pendent parameter is determined by index K4 of PA(K4), where K4=l 
for nrp, K4=2 for T^, and K4=4 for T^, The program selects the 
dependent parameter in accordance with the profile, either PROF=l 
or PROF=2o 

FORTRAN symbol 
K4 

NPA 
PA 1(10) 

Output 

An example in appendix J shows the input data in the same format 
as punched on cards for four wavenumbers and two gases, H 2 O and N 2 
The output has four parts: 


Description 

Index of independent variable temperature parameter 
PA(K4) 

Number of values of PAl to be used, up to ten. 
Values that may be assigned to PA(K4). 
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1. The NAMELIST aad READ data 

2. Temperature and pressure profiles / 

3. At each solution found by iteration, the four variable profile param- 
eters, wavenumbers, number of iterations, final r ms residual 

4. Tabulations of dependent parameter against independent parameter 
for each solution in (2), with the dependent parameter modified for 
plotting by equation (27) with the constant ^ 

Examples of plotted tabulation are figures 3 and 4. 

DIAGNOSTICS 

A warning is printed each time: 

1. The pressure of all absorbing gases exceeds the static pressure 
at any point on the pressure profile. 

2. Ehiring iteration when a variable parameter reaches the assigned 
upper or lower limit. 

APPLICATION OF THE PROGRAMS TO AN EXPERIMENT 

Three steps are used to determine whether this method will be suf- 
ficiently accurate for an application. These steps will be summarized 
in this section, and an example of an application will be given. The steps 
are: 

1. Estimates of the temperature and pressure profiles are made, 
and the profile equations are selected. 

2. Some wavenumbers for measurem.ent are selected, and then com- 
pared using graphs of the output of the WAVES program. This was 
demonstrated in figures 3 and 4. 

3. Accuracy of the dependent profile parameters is determined for 
the best pairs of wavenumbers selected from graphs in step (2). This 
was demonstrated in figures 5 and 6. Tests with the PROFIL program 
in MODE -3 may be run with a perturbation of an independent profile 
parameter. The resulting change of all dependent profile parameters 
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may be judged acceptable or excessive. This is demonstrated in ap- 
pendix I. 

based on the results of steps (2) and (3), it can then be determined 
il there exists a pair of wavenumbers that lead an acceptably low b^vel 
of random error in the computed parameters. The experimental meas- 
urements are then made at the wavenumbers chosen and the results 
serve as inputs to the PROFIL program in MODE =2 for determining the 
temperature and pressure profiles. 

An example of the PROFIL program is given to illustrate the limita- 
\ion of the assumed functional form of the temperature profile as de- 
scribed by two variable parameters n An assumed profile with a high 
temperature core xd shown in figure 9. This profile cannot be matched 
exactly by equations (la) or (lb). The pressure profile is the same as 
used in other examples; hydrogen-air combustion at 1 atm, parameters 
p,„ = 0. 30 atm, p„ = 0. 35 atm, and n„ = 0. 50. 

W C ^ p 

By following steps (1) to (3) given at the beginning of this section U 

was found that wavenumbers of 3100 cm"^ and 3500 cm"^ were suitable 

for use. The radiance and transmittance of the assumed profile were 

calculated at these wavenumbers to represent error-free measurements. 

With these measurements as input data, the PROFIL program in MODE “2 

found the profiles with equations (la) and (lb) that are plotted in figure 9. 

The profiles in figure 9(a) and (c) use an independently determined wall 

temperature T„, = 1100 K that was made 100 K lower than the assumed 

profile T„, = 1200 K. This represents an error of T The difference 
w w 

between the assumed and computed profiles is an error of approximation 
due to limitations in the ability of the profile functions to assume an ar- 
bitrary shape. This may be compared with the error band tliat is due to 
an assumed radiance measurement error of one percent and a tabulated 
absorption coefficient error of five percent. 

The PROFIL program data cards and output for figure 9(a) are shown 
in appendix K. The output format is the same as shown and described for 
MODE=3, appendix I, without the first namelist. 
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CONCLUDING REMARKS 


/ 

All of the examples presented in this report were for hydrogen- air 
combustion at a pressure of 1 atmosphere with a path length of 100 cm. 

These reactants produce H2O as the absorbing gas, IL was shcwii in 
table in that the random error in the computed profiles decreases with 
increasing path length. It can also be shown that the random error in 
the profiles will decrease with increasing pressure or by using a gas 
with greater absorptance such as COg. However, in applications where 
the path length, gas pressure, and gas type are all fixed, a user of this 
computer program has only the choice of the wavelengths of measure- 
ment at his disposal in order to attempt to bring the computed profile 
errors within acceptable bounds. 
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APPENDIX A 
SY.MBOLS 


A 

Ai to Ag 

a 

B 

b 

d 

k 

L 


Hrji 



m 


N* 



P 


Ps 

q 

R 


R (subscripted) 
6R (subscripted) 
R (subscripted) 
6R (subscripted) 


matrix of first derivatives 
coefficients in eqs. (16) and (17) 
element of A 

matrix of second derivatives 
element of B 

average spectral line spacing 

gas absorption coefficient, cm“^ 

total path length of profile 

exponent in equation for temperature profile 

exponent in equation for pressure profile 

actual gas radiance 

temperature- gradient, eq. (2) 

radiance by Planck function 

apparent radiance defined by eq„ (24) 

partial pressure of absorbing gas 

static pressure 

number of absorbing gases 

total residual, rms value; also column matrix of 
residuals 

residual, eq. (6) based on and 
residual errors, eqs. (20a) and (20b) 
residual, eqs. (23), based on 
residual error, eq. (25) 
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T temperature 

temperature parameter in temperature profile, eq. (2), fig. 2(b) 
u dimensionless optical depth 

V detector output voltage 

w independently varied parameter 

Xj boundary layer thickness 

X distance along profile or optical path 

y dimensionless oath- length ratio, eq. (1) or variable parameter, 

eq. (8), (Bl) 

Y graph- plotting ordinate, eq. (27) 

Z dependent parameter 

55 dimensionless path- length ratio, eq. (2) 

absorption factor of gas 
Q!* spectral line -broadening coefficient 
/3 spectral- line overlap factor, band model parameter 
y average spectral line vndth 

7] interpolation parameter 

X wavelength 

p gas density 

transmittance 

^ graph- plotting constant, eq. (27) 

Subscripts: 

c centerline of profile 

s standard conditions 

w ends of profile (wall) 
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p pressure 

T temperature 

N radiance 

e effective value 

L total path length 

a absorptance 

X wavelength 

0 initial conditions 
Superscripts: 

c computed 

m measured 

Running indices; 

1 path length; also wavelength in eq. ( 6 ) 

3 absorbing gases 

k nonabsorbing gases; also dependent parameters in matrix solution 
L intermediate index for path length 
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APPENDIX B 

PROFILE PARAMETER ERROR ESTIMATE 

In order to use matrix notation and FORTRAN coding, the four vari- 
able parameter symbols T^, T^, p^, and n,p are replaced by one sub- 
scripted symbol yj. Only three of these are treated as dependent vari- 
ables at any one time. 

Due to random errors in the measured and computed terms that 
form the residuals and R^, there will also be a random error 
6y^ associated with each parameter y^. The random errors in the re- 
siduals are denoted by (see eqs. (20) and (22)). 

When the number of equations is equal to the number m of unknowns, 
the errors 5y./y. are a solution of m equations in m unknowns, 

J J 



or, in matrix notation, 

A^^j:.6R ^B2) 

The elements of A are partial derivatives in equation (Bl) that are founr. 
by numerical differentiation in the computer program. 

When the number n of equations is greater than the number m of 
unknowns, the least squares method is applied to give the equation 

= (B3) 

(■p 

where A is the transpose of A. This can be written as 


/ 
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(B4) 


which IS equivalent to the set of n equations, and 


m unknowns, 






6R. = — i 

1 


(B5) 


Each term in equation (B5) isolates the contributions of to the 
total error of each variable parameter. Since the sign of 6 r^ is 
equally likely to be plus or minus, the rms total for all n terms gives 
the estimated error for each j parameter. 


Damped Least Squares Method 

Where an iterative method of solution is used to deduce the profile 

parameters from experimental measurements, the profile parameter 

fractional increments a-*'e ^Yj/Yj for an iterative step. The increment 

size is determined by the residuals R-^ which are to be minimized. 

These residuals are: Rj^, equation (23), at each wavelength X^; R^, 

equation (7), at wavelength X^^. The least squares method, equation (B3), 

minimizes R by reducing ATr to zero. Near the solution where 
T 

-A R becomes zero, dammng is required to prevent oscillation of the 
variable parameters about thr solution (refs. 17 and 18). When a damp- 
ing factor is selected in accordance with equation (11) of reference 19, 
equation (B3) becomes 

|a'^A + (p/q)Qj ^ = -a'^R (B6) 

The elements of matrix A are first derivatives of the residual, given 
by 



yj(3Ri/3yj) 


(B7) 
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where R. is by definition a fractional value. The transpose of A is 

T 1 T' 

A . The damping increases the principal diagonal of A A by a factor 

1 + (p/q) where Q is a diagonal matrix whose elements are the 

principal-diagonal elements of A A. The number q is the length of 

T 

the principal-diagonal vector of A A whose elements are q^ ^ so that 



The number p is the length of a vector given by the matrix product 
T 

B R, whose elements are p^ so tliat 



The elements of B are second derivatives of the residual given by 

by = y?02R./3y2) (BIO) 

Finite differences are used to find a^j and b^j. First treat the profile 
parameter y^^ as follows for all wavelengths: 

(a) Calculate current values of Nj^ and r-^ 

(b) Increase y^^ by multiplying it by 1. 01 

(c) Calculate new values ^ and ^ 

(d) Decrease y^ by a factor 0. 99/1. 01 

(e) Calculate new values Nj^ and 

(f) Restore y^ to its starting value by multiplying it by 1/0. 99 
The first derivatives for the ith wavelength, jth parameter are 
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^r<X,b-X,a)/C0-02(l-”NL)] 

^j=(“%b-'%a)/C0-02(l-“rL)] 


the sec.ond derivatives are 


>’ij = -<X,a ^ - 2X>/C0' 01^(1 - “NlO 


blj = -('’’-L.a - '^%b - 01^(1 ' ”-l)] 


(BID 


(B12) 


Next, repeat the procedure (steps (b) through (f)) for the succeeding pro- 
file parameters y,. 

J 
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APPENDIX C 


DETERMINATION OF ? IN EQUATION (27) 


The WAVES program involves the determination of how wavelength 
A affects the relation between a dependent temperature profile param- 
eter Z and an independent temperature-profile parameter w. At the 
initial value w^, the value of Ztw^) is invariant with wavelength. How- 
ever, Z(w, A) changes slightly wi.th wavelength when w w^. The 
temperature-profile determination depends on these small changes. A 
graph of Z against w on an ordinary sheet of graph paper would show 
curves of considerable slope, so that differences between curves for 
various values of A would not be readily apparent. To exa.ggerate 
these differences, the ordinate is changed from Z to Z + C(w - w^). 
The value of ^ is determined as follows. 

Let the range of the independent parameter be 


w, 


n 


w s w, 


m 


Among the values of Z(w) for all the wavelengths which are tried in 

the computation, there will be a maximum value Z^^^(w„) and a mini- 

iii^x n 

mum value Similarly, there will be a maximum value 

^max^^m^ and a minimum value Zj^in^^m^* value of ^ is taken 
as 


-Z (w)+Z.(w) Z (w)+Z.(w) 
1 max' n-^ min' n' . max' m'^ min' 

— — + — 


(Cl) 


The total span of ordinate values is thereby reduced severaUold. 
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APPENDIX D 
PROFIL PROGRAM 


Program Symbols 


The following FORTRAN symbol list includes symbols in the sub- 
routines of appendixes E and F. It does not include symbols defined 
for NAMELIST and READ data in the section PROFIL PROGRAM USE. 


FORTRAN symbol 
A(5,5) 

B(4) 

BB(8,3,41) 

BE(8,3,41) 

BO(8,3,21) 

CP(6, 41) 


Cl 

DAMP 

DEL(4) 

DPA(4) 

DPAPA(4) 


Dc cription 

Condensed matrix of first derivatives 

Condensed matrix of second derivatives 

Denormalized and broadened BO 

BB weighted by KK 

Interpolated local value of BTAB 

Constant of pressure profiles; ratio of each 
absorbing gas pressure to first absorbing 
gas pressure; ratio of each nonabsorbing 
gas pressure to the total pressure of non- 
absorbing gases 

Constant of pressure profile of first absorb- 
ing gas, PW(1)/PC(1) 

Damping constant for diagonal of least 
squares matrix 

Fractional change of variable parameter 
for derivatives (in data statement) 

Variable parameter increment for iteration 
step 

Fractional value of DPA 
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FORTRAN symbol 
DPV(3) 

i)TAU(8, 41) 

DUM 

DUMl 

DUM2 

DX(41) 

ERR(IO) 

ETA(8,3, 41) 

G(10,6) 

GSAV(8, 4) 

GV(9,5) 

H(9,5) 

HV(9, 5) 

lA 

IB 

ITERl 

12 

JMIR 

J2 

J3 

KO(8,3,21) 

KEY 


Description 

DPA in condensed matrix 

Increment in transmittance 

Dummy to prevent overflow or underflow 

Dummy to prevent overflow or underflow 

Dummy to prevent overflow or underflow 

Distance between stations X along path 

Residual caused by errors ERRV and ERRK 

Interpolation parameter 

First derivatives for least squares matrix 

Dummy variable equal to G 

Dummy variable equal to G in condensed 
matrix 

Second derivatives for damping least squares 
matrix 

Dummy variable equal to H 
Index parameter equal to I 
Index parameter equal to II 
Iteration counter 

Wavenumber index parameter I+l 
Index for reflection of profile about center 
Path index parameter (Jl+l)/2 
Path index parameter Jl+1 
Interpolated local value of KTAB 
Selector for variable parameters 
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FORTRAN symbol 
KEY1(8) 

KEY2(8) 

KEY3(8) 

IOC(8,3,41) 

KV 

KV2 

K1 

K2 

K3 

LX 

L3 

MAX 

NEWLIM 

NPl 

NTEST 

P(6,41) 

PAERR(4, 10) 
PA1(4) 


Description 
Index parameter (data statement) 

Index parameter (data statement) 

Index parameter (data statement) 

Denormalized KO 

Index of condensed matrix 

Index parameter KV+2 

Index parameter for variable parameters 

Index parameter for variable parameters 

Index parameter for variable parameters 

Index parameter for error array 

Index parameter for gas type 

Largest absolute value of all DPAPA for an 
iteration step 

Controls number of tests for upper and low 
parameter limits 

Equals pressure profile equation exponent 
NP(1) 

Diagnostic indicator of absorbing gas pres- 
sure being greater than static 

Pressure of each gas at each station along 
path 

Error of variable parameter caused by a 
residual ERR 

Tentative value of new PA for an iterative 
step 

Sum of pressure of all absorbing gases along 
path 


PSUM 
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FORTRAN symbol 
R(8,41) 

RA(8) 

RAB(8) 

RB{8) 

Rl(8) 

R2(8) 

STEP 

T(41) 

TAU(8,4, 41) 
U(8,3, 41) 
WAVOUT(IO) 
X(41) 

Y(21) 


Description 

Gas radiance 

Equals CRL during numerical differentiation 

Equals CRL before numerical differentiation 

Equals CRL during numerical differentiation 

Constant of Planck equation 

Constant of Planck equation 

Sum of fractional changes to variable param- 
eters for an iteration step 

Temperature at each station along path 

Transmittance at each station along path 

Optical depth at each station along path 

Wavenumbers in sequence for printout 

Actual path lengths Z* PATH/2 

Dummy equal to Z, the normalize:] path 
length 
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PR0FIX<E Program liBtdng 


Lbl^l 

1* 

..OTOl 

?* 

IbtOl 

3* 

UblOl 

AA 

bcni 

5» 

i oni 

b* 

bbl^O 

7A 

LblOM 

AA 

bOlO'j 

9A 

bOlOb 

lOA 

UblOb 

11* 

Obi Oil 

12A 

uc 1 0 

13* 

UOlOb 

lA* 

b01D7 

ISA 

ccir>7 

IbA 

00107 

1 7* 

Lono 

1 BA 

bUlIO 

19A 

uOllT 

20A 

bUl 2vi 

21* 

00123 

22* 

Lbl?b 

23A 

bbl3b 

2A* 

bOlM 1 

25* 

..bl'b 

2 b* 

OblTl 

27A 

t0174 

2B* 

lb 17b 

29* 

UOOOu 

30* 

bU203 

31* 

CD20b 

32* 

00211 

33* 

bU21b 

3AA 

UD217 

35* 

00221 

36* 


y.2n.t3. 

X BTAb(e,J,B),bOlf.,3,2n,l'Pie.3;2l> .^tia,3,4n,tT4|B,3,MU, 

X U(B,3,«niTAUtP.«i‘lll,NPl»IAilBtNTEST,TX 

PCAL XX ,XT4B,KC,iNPl 

DHE5riprnpIpMA>.PAP<p^0P*'-''**<‘-^i'fi>''<^"2(8),KCY3.e),Daix^ 

X HT»UU(RI.CIlD,Pl,H19,5UOV(9,51,HVt9.B1.A(5.5 .PAEPPtS 

X , 3 «Sfbl,W«V£;i 6 )iZ (21 >tr»PllJ>»XAVOUT)nt ,CSAV(B,H,PW(b), .PUI, 
X PC(b) 

9AXE ulSI/NAHl/POaE!pBOF,IBt;An|PATH,n,WAVE,HPL)MTAUU. 

X EPPvtrPXX^w;K,NT;TU3,:.B.UPL:«,POTU,H,L3,l.2,PS,PC.MP.P-,NB. 

X TKTAi',Ni3,TPTAB, JliZ,ITERfTX ^ , , , , . „rvxyi 1 2.MA1.U/ 

•)ATA XEY1/4*1,2,2.3,D/iKLY2/3,2,3,1i3,2,3,3/,KEV3/1 ,1 . 

’C UEAO:a 

Pi.APt‘’.'<A‘«I ,EN0:9?9) 

IMIsrAD.FOiOIGO TO b7 
9EADtb,BC,t(GAbU l,L = l,bl 

PEAn(5,nn ICALPP(LfLL).U:i.M ),L-1,U) _ 

PEADCi.PlMIIKTAtin.LtNI.NPliP'tlPl'IJ '*-:*'^' 

PEAO(5.Sl)n(PTABniLtSltN:ltPI»lPl?ll'.L-l(LI) 


00 bb L = 1 ,H 
alpH(L.L>p1. 

00 bb I = lf I 1 
00 bH V=l,NK 
frft [F(KTAB( I ,L iNI.LT . 

00 bb R=ltNB 
bb IF(PTABa,L,NI.LT. 
b7 CORTIWUr 
XRITEtb,90l 
XBITE(b,NAHl) 


.H -Obi KIAO ( I ,L »N> F-lt-Ob 
,lE-0blBTA"(I»L|NI = .U-3b 


uDon jO 
LOuniiD 
UOpPljO 
bOunuD 
bCOOjO 
tiOU^UU 
amiPbi 
liOJPCil 

ononoi 

QDUOUl 
UQOOOl 
uobpoi 
liDjroi 
□OOPUl 
uOtinbi 
LPJPSl 
bOJ^J 1 
bPOPbZ 
b0bPb2 

bnopoz 
b'>bOb2 
bPb''U7 
tnb'’n 
uDj^MS 
uCLf’b7 
bOull 3 
bPulHb 
0''- lAB 
bOulbCI 
uHcibO 
L'^blbb 
bOJ IbJ 
bOulbl 
tObPuO 

m0u'’00 

bnj?LA 


ORIGINAL PAGE IS 
OF POOR QUALITYI 
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I-LOOm 



If 1 JwrA'^.rO.OlG'’ TO 23? 


tp* 


rfMTOlt.bl) (2A5 IL 1 ,L:l ,t 1 

. 12M 



2MtL(b ,«7 1 t ( ALPHIL ,LU tLL =1,41,2 = 1,21) 

u.:‘i2 



KPI TE (6,9“ 1 ( ( (KTA(>( 1,L,N1 ,N=1 ,B) ,1=1 , 111 ,2 = 1 ,2 1) 

;,i 2«t, 

U 1 « 


4HTt(6,b9)(((bTArtn,2,M,N=l,b),I = l,lll,2 = l,211 

t 025 l 


c 

IMTIA2 CONDITIONS 

io:7^ 

4 3* 


VTE5t=G 

1 •'i 

44* 

500 

IA = I 


4C,* 


I [-=11 

! 02’’b 

4b* 


0CNE = '’ 

tl.27b 

47* 


I TEX 1 =n 

uo;"’? 

4P* 


0»4p:o. 

v030a 

49* 


i 2 =n*i 

>.:,>ai 

?n* 


J3=jm 

1.2202 

51* 


J? = ( J1*I 1/2 

.'utOi 

P?* 


. 3 : 21 *1 

; 03Dm 

53* 


C!:'>4| 11/DCIl ) 

i.D30b 

54* 


PA( 1 ) =NT 

litjJOt 

55* 


PA(tl=TC 

UC;3'17 

Gfi* 


PA ( 3) =PC( 1 ) 

233IJ 

57* 


PA( 4 1 =TW 

(.2 311 

5P* 


DC 64 K=I ,4 

02314 

59* 


DPAPA(K1=0. 

2231b 

bC* 


3PA(A)=D. 

uC2tl, 

61* 


CTAUtlK )=". 

12317 

6?* 


:tau2(k*i )=o. 

1.332o 

fc?* 


DO 64 1=1,10 

(.(,'32 3 

64* 


/AV3uT(n:ri. 

o2324 

65* 


EB4(I)='7. 

u037b 

tb* 


5(1,6)='’. 

UJ37t 

67* 


3(1 ,5)=D. 

2L.327 

6P« 

2 4 

»AEPP(6 , I 1 = 2'. 

o232i 

69* 


DO SO 1=1,11 

u23''l. 

7C* 


IF(MDDF.tC.llCO TO 65 

Ijj23 7 

71* 


WAVOUTl t 1 :W44E (I 1 

wk. 3 4 u 

7?* 


WAVOdT(K') = HAVL(l) 

(■.214 1 

7?* 

65 

41(11= .119D9C-11 *2AVEin*93 

22342 

74* 

51 

421 1 1 :1.4 38B»WAVt n 1 

CC344 

75* 


DC SI J=1 , J2 

23 34 7 

76* 


JMIP:J3-J 

i>2'3'2 

77* 


XlJ> = PATM*2(JI/.7. 

u23' 1 

79* 


V(J1=2(J1 

(.0322 

79* 

51 

K ( J«IB) =P4TH-X ( J) 

(10354 

BP* 


DC 55 J=2,JI 

23357 

61* 

59 

DX(JI=X( JI-*( J-1) 

2L357 

8 ?* 

C 

PXESSU4r OPOFILF IN1HA2 PBOFTOfS 

23321 

63* 


NFUNP) 1 ) 

2U3(>2 

84* 


DO 52 3=1, J2 

20365 

85* 


DO S3 2=1,22 

2D373 

66* 


DUH1;APS( °C (2 )-P>, (LI 1 

22371 

R7* 


PI2,J)=0C(2) 

20372 

88* 


IF(CJ«l.tO.C..OO.J.E0.J?)50 to S3 

22374 

89* 


DU4 = l.-’'( J| 

00375 

9n* 


DU4 2 = A2D;. IDUM1)*AL02(DJ«I/NP|L) 

U(j37b 

91* 


3F(0U“2.27, -1-2.120 TO S3 

2C40J 

92* 

5 6 

PlL,Jl:°CfLl-(Pr<LI-Pw(2) l*0u“*9( 1 , /\T (L) » 

22401 

93* 

53 

C ON T Ivor 


uOi-’ir, 

L'luZiJ 

k."u?£,2 

t,n;j’’7f. 

uTu’Th 

OHU722 

L,nj’2« 

bHu’t? 
LTii’JC 
Lnii’Jl 
COL’32 
Qnj’iu 
Lnt’jfr 
0'’u''HC 
u'ltj3Hl 
C'^iJ ■'<(6 

uriJ3:.2 

Lnu3i,K 

u''i-117 

t,fiD«17 

Lfu'in 

cr’cu 20 

uCu*j2?' 
unuijti 
"uUOS 
uOu 2(> 

0'’u‘<3: 

wOJUl,? 

LOuI-mO 
ti'iuxM 3 
UJui-Mb 
t-r’-HtO 
Lr;iiti4 
2 

fc*n j M 1 6 
Lf';,M72 
00J><7M 
LnO''L‘( 

U^u'U^ 

LO„'m 

:-0^'.u7 
uOj‘ lb 
[jOl'- 2f- 
LOi.tJt 
OOi'Ml 
U'’u'>t2 

L0u‘’1.3 

L"u'tS 

i.0L,''7C 

LO^MC 
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.,q?3 

bq* 



iriLi«70,nriO to 55 

..tiqP3 

9** 

c 


»B£SSi'qr PooFitr bmio 7*ctohs cpil.ji 

i jqjb 

9t» 



30 £q L=2,Ll 

f blq 

97* 


r q 

OML 1 J):P(L ,JI/Pt 1 1 J) 

CM 12 

9A* 


7 7, 

IF (1,2.70.1.1 )(,c TO 52 

q^.qiM 

SO* 



PUh:J. 

.,1 q I b 

1?0» 



30 5b L71.7,L2 


ni* 


7t 

Duq:aiJM*P(L,Jl 

.it,q:i 

J32« 



DO 57 L9L3,L2 

ubqrb 

1 n''* 


S7 

CPIL,J|:?IL.J)/njM 

. t-qr? 

Ijq* 


f 2 

CflNTIMUF 


n'» 



IFlMJOF.tC.DICO TO 731 

i.:q33 

nb* 



CALL 7fqp 

q( q3q 

1P7* 



C»LL 70*95 

qiq’b 

1 3P* 



DC 63 1=1, n 

i b’qqj 

no* 



qPLlt nCRLU ) 

..uqqi 

1 in* 


f 3 

OUULtl 1 bOTAJLI 1 1 

w^qqi 

1 1 1* 



! 71 poor .60. 1 K'O TO 701 

i.k,qqb 

112* 



qODC 

; Lqqb 

1 1 3* 



50 TO 505 

bStq? 

I lb* 


721 

2M16l6,9niPMII,I:l,qi,Pl.(lt|NPai 

,uqq7 

11'* 

c 


36319 ITL3*TI0N loop 

vCq'.q 

lib* 


7rr 

C*lL T6PP 

Lqqbb 

117* 



f*LL tOA'ib 

„iqbq 

IIP* 



I7| 9ir5T.ro. 1 |wPlTF 16, 731 

. i,qhi 

110* 



9T7SI :o 

i-bqfci 

12 3*' 

c 


7190 qrSTDUALS, TLSI 709 SOLUTION ON MAXTHUH 1T6PAFI09 NUMPf P 

„oq' 2 

121* 



9651 J:n. 

bqq3 

12 2* 



r AJAit:CTAUL ( 1 1 

bijqbq 

12^* 



3C 1 191,11 

tqh7 

l?q* 



9*3(11 = 791 { 1 1 

:q7j 

12b* 



5(1 ,5I=C9UI I /M9L( Il-l. 

q7i 

1 2b * 


1 

9F515 = «r5I0<>0(l, 51**2 

<..q7j 

1 2 7* 



If (“J"6.tO. IM,0 TO q5b 

uq7a 

13P* 



5(I2,'l=(‘'TAJLm-CT AJLtl 1 )/( 1 .-MTAUL (11 1 

oi q?:, 

129* 



96SIJ = '.09T(B651"*7, (I2,5I**21 

„q3^ 

1 1** 



5T6P:’'. 


m* 



3C 2 J=1 , 7,2 

‘ 3 -i 

172* 



9I7p:5irP*AH5(LPAPAlJIM.l 

w'^q 

177* 



SUP = 7Tr>,*PS (5PAF A ( J*1 1 1 

•bbPb 

1 3q* 



I f ( ( 1 T69 1 .OT.C 1 .AND. ( ST7P.LT. . 5031 .09. 17L9 1 .60 .11 79 .09 .P65 P.LT 

. C 'Tb 

17q* 



X .3331115796=1 

‘ 1« 

1 3b» 



5911616, 57M0P*(I1, 1 = 1, ql, IT L7 1,9rSI0,5T6P,(PAm, 1 = 1, q), 

i .■'Iq 

1 3 7* 



X (31 I ,51 , l:l,SI , (WAVOUT (I 1 ,1 = 1 ,51 

-wb’l 

I 7 P* 



I T7Rl:ntPl*l 

' Jb?! 

1 39* 

c 


3f 91VAT7VF S 65(311, 3121 

.b‘.?2 

m* 


4 ?( 

<f T = t 

.* ** 

iq 1* 



3l:K671Urv)-l*np07 

L..'*2q 

lb2* 



<2= KrY2(«6Vl-l*P907 

vbTb 

1 q 3* 



<3 = '<6V3HrYl 

i b2b 

1 qq* 



31 .( 1 1 9.D1/PA ( 1 1 

,...b?7 

Ibb* 



DC 3 5:< 1 ,k2,93 

q ^ t. 

1 q f.* 



p A( K 1 rP A 1 K 1* ( 1 . *06 L (9 1 1 

> 

1U7* 



OA.L TFHP 

i bb3q 

1 q «* 



CALL TPA9' 

b‘. 3 -.I 

1 q 9* 



30 q I = 1 , r 1 

.,{ *qq 

1 7 P* 


d 

9 A( 1 1 :C9w 1 T 1 


LPoMO 

uru^ J2 

LPu* n2 
C-"-f £-1 
u<'0'-bl 
u'’uf'bC 
w.lUf'bO 
U'lU'ibJ 
wOj'bS 
0'’Ch72 
LnL672 
bnjb73 
uCG h 7 b 
u'7U7Lj I 
L''u’L3 
bCti 7 t;b 
bTJ‘'36 
L'7b7tG 
b0u7,;l 
j'’j723 
U'1J7J3 
1 3 3 
b'lj’bM 
cnu’35 

U"l<7qb 

tnii’qb 

bbu7M6 

L'lu’bO 
uP J’bfi 
t."b'’7S 
I ni'’bl 

banjb 

oPl"lU 
tniPiD 
bP 1 q 
[ n 1 "qq 
bPl''7l 
uPl"'71 
I." 1 " 7 b 
i.Pl"7b 
l.Pllq3 
L,n U7 

i-''nn 

c'’l''3q 

wP^ql 

b^i^q', 

w'’»pq7 

u"l ’bM 


ORIGINAL PAGE IS 
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. u'ib.' 

1?1» 


rAJ«:ri«oL 1 1 1 

,.l:' b •. 

I'r* 


PAIK)-P*lK|*M.-D('UA||/(l.»LnLMI) 

..u'jbb 

l'J» 


CAww ’IMS 

, u‘ ML 

I'M* 


CA.L T9AVS 

lU'-Mt- 



30 5 1:1,11 


I'i7 * 

<, 

TMiiTCs.m 

uu'l'L 

l';7* 


lAJprCTAuLM ) 

U'>' M 

I'.ft* 


»AI A) :daIb l/( 1,-UEL (K ) I 

lU'j'L 



36 6 i:i,n 

'.u66l 

itn* 


3 1 1 ,< 1 ; ts A 1 n-961 1 ))/ 1 *nu iH i*mi,l i r n 

u'lH 

IM* 


:,'Avii,i<»-on,Ai 

1 'i 6 fc 

i6r» 

r 

M(I,AI: (5Alll'*9ii(tl-r,*M*n(lM/(ri£U>il**?*''»Lll)) 

t 'j6b 

163* 


j( ir,M 1 ■( 1 AUl-lA J* J /( ?,*3f lIA 1 *(1 .-MT AUL( 1 ) 1 1 

uU!)6L 

16.M* 


MU r ,«>:l 'AUA7I A J3-?,*1A JAM) / IU( L Ik ) **3*1 1 ,-Hl AJL ( 1 1 n 

i.U6 6b 

It '♦ 


Ct/MI Ibun 

i U 'J 7 u 

166» 


CA.;. 7LM3 

1 - 7 1 

167* 


CAUL tSAS'. 

.U'7u 

160* 


S£ WU1«71 

,,ub ’ :■ 

IfeM* 


36 J2 K:1 ,M 

1 3' 7b 

173* 


■’AKUrPAUI 

u67 / 

171* 

7 *7 

-IHlKCn 


17,’* 

C tn* 1537:1, LCUUtM I4LL7 6 hI6T IlRIVATIVIO ADOVL I:.’ TO 

u' 11 

17 3* 


•jc Mil «iT.r,ii 

. . (."b 

17lj* 


If (MLnr.sr.iibc in ‘.nn 

: <i'’j 

17',* 


ir:j 

. 363 7 

1 76 * 


aAV.IJKI l-wAVf 111 

. l - 

177* 


WA9Lj7(’)rwAVE (“1 1 

u6 1 1 

1 7fl* 


AA\/5j7(7):i,nyt , j) 

i blu 

1 79* 


[ 3 M’i K7Kl,Ki,K3 

u61‘. 

i»n* 


jl3,U:n| IH1,K) 

L' lb 

ini* 

A ' t 

j 1 r I A 1 :r.LAv (M 1 ,B 1 

• ta W 

1 «.’* 

f, -T- r 

ccvii'lur 

f 

1 ' 7* 

C 

CHAvior niLlMAllVL VA9J49I.I.5 in CONOEI'-t A99AV K:Kl,Kr,K.3 in K 


IB** 


A V:3 


1 1. f, ♦ 


30 7 k:k1,k2,kJ 


106* 


AV:av*1 

t ZhPu 

1 07* 


30 7 l:l,I2 

. *» M 

1 1- U* 


3VI I ,KVI :f (I, M 

u-b’.. 

I 0 9* 

7 

AVI I, A V) --HU,*) 


1M3* 

C 

.lALT SCJATtO MATmJX AT--A [w(o6| 

■ A 

ni* 


30 0 i:i ,KV 

.'PfUj 

1 9?* 


3C B k:1 ,«V 

.. :»-M r 

197* 


A II , A U P . 

^ b *4 M 

1 9M* 


30 f j:l ,12 

i,/ (4 / 

196* 

fi 

Air,K)AA(I,KH6VIJ,KI*GVUl,ll 

...t bi 

196* 


rnM0nE.c,Cull60 10 M53 

'.3 

]97* 

c 

9L6I0UAL VECIOM C,T--M £C(°6I AND nA“PINr. VOCIOH 3 f'lPbl 

..-.b'L 

19S* 


<V2:AV*i 

:,7'b 

179* 


16 9 K:1 , BV 

. : 'f 1 

.’13* 


A(A,AV?):ii. 

. 1 ' b J 

rri* 


0IA|:i. 

r<‘6 ; 

?3 ?* 


30 9 j:l , T2 

u ' 6 b 

.”7* 


4(<,AV2):AI‘',Kl/’)-n,V(J,B)*r.lj,L) 

„bb 7 

r3M* 


3(A|:'MA)*hvIJ,B1*GIJ,5I 

L67^ 

?C 6 * 


damp:’. 

l67 / 

,'no* 


IF (i 2 .ro.Kv.aB .OGNC .cc. nr.D m ip 

-67^ 

33 7* 

c 

VCMwALI^rn damping vrcinp P/„ iOlBtl, P9I AVt UAMnc’ maimj* 


C.Ti-’iib 

I.'n.'’b7 

u"V-lt 
uni’/b 
LTl »ijC 
I "1 7u? 
u:'l »l<! 
L"i Mr 

(i''l''17 
I ”! M l 
uiiniJ 
L''l M2 
wtl 1 ^b5 
bill lb', 
uni '67 
onnii 
uni<u7 
Lniq^'l 
031<130 
C.niMiC 
f.niqjl 
uri‘((>i 
unii(63 
i.ni<tb5 
bt lHfa7 

unio? 1 

miriL, 

uni' 1C 

-'U> 1'. 

uni' lu 

uni' 1 ' 

U"1'U3 
l.•'l'b^ 
uni 'Lb 
uni' t 3 
uCi'm 
U^l'm 
uni' 76 
^ni> 1^ 

U n 1 6 1 
n 1 M (. 

U"! ' Ifc 
uni'37 
uCl' 37 
u^l'MC 
uMbtL 

u"!' jb 

uni^'b 
u"!' br 
I '’I'br 

u'l'b' 
-•’iTj: 
uni M7 
Lni 7^3 


ORIGINAL 
OF POOR QUALITY 



50 


V w 7 ‘J 

!V>* 



0U3. 


U0J712 


ro9* 



32£3. 
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APPENDIX E 

TEMP SUBROUTINE LISTING 


1,0101 

1 * 


SUBROUTTNr TEMP 

.,6103 

2 * 


C 04 M 04 0405 , TWb, 0 1 , J 1 , 32 , J 3 , C 1 , PSU M , JB , R 1 ( B ) , S 2 ( B ) , X( 4 U , 

00133 

34 


X F'fl ( 4 ),CT 4 UL 16 J,C 4 L( 8 ) ,T |4 1) , 0 X (41 ) ,NK,NB,PS,pATH,LZ, V (21 >,L 3 

00103 

44 


X LT>UlB, 4 n, 4 tb, 4 l),CP( 6 , 4 W,P(b, 41 > , ALPHl 3,6 » , 

.0103 

54 


X TXTAOI 8 ) ,KTA 3 ( 6 , 3 , 8 ) , KOIB , 3 , 21 ) ,KKta , 3 , 41 » ,T 3 T 46 ( B) , 

.10103 

64 


X prA 3 ( 6 , 3 , 3 ),HQ(B, 3 , 211 , 9 B( 8 , 3 , 41 ),PE( 8 , 3 , 4 l),ETA(B, 3 , 41 ), 

60103 

74 


X UIF,I, 41 ),TAU 1 P, 4 , 41 I,NP 1 ,IA, 1 B,NTEST,TX 

00104 

64 


4 EAL KK,KTAB,KU,NP 1 

00106 

94 


I 4 TEDER PROP 

6 313 5 

134 

c 

TE 4 P A 40 P 6 r 4 CIPAL P 6 ESSJPE P 90 FILC EO(I , 2 , 3 ) 

60106 

114 


PUl :PA( 3 ) 4 C 1 

00107 

174 


DU 419 AB 5 (PA( 2 )-PA( 4 )) 

OCI 16 

134 


0 U 42 rARS(PA( 3 )-PWl) 

6011 1 

144 


00 1 3 = 1,32 

03114 

154 


34 IP= 33-3 

ODl 15 

lb» 


T) 3 ):PA( 2 ) 

00116 

174 


DU 4 - 1 -Y 13 ) 

301 1 7 

164 


IF( DUMl.EO.O, . 0 P. 3 .EQ. 37 IG 0 TO 11 

C 0121 

194 


Du 53 :AL 0 G(CH)Ml)+AL 03 (DUH)/PA(n 

JC 172 

204 


IF(DU 43 .uT.- 8 C. )G 0 to 2 

60124 

214 


11 T( 3 )=PA( 2 )-(PA( 2 l-PAt 4 ) ) 40 UM 4 * ( 1 ./PA) 1 ) ) 

3 C 1 Z 5 

724 


7 IF(y) 3 I.GT. . 51 ) T( 3 ):T) 3 )«TX*lSlNt 3 . 141 b*|Y( 3 )-. 5)))**2 

0 C 1 Z 7 

Z 34 


T( 3 MIP) PT( 3 I 

GDI 30 

244 


P( 1 , 3 > = PA 1 3 ) 

00131 

754 


IF) 03 “ 2 .EO.O. .OP. 3 .E 0 . 32 )GO TO 1 

uCl 33 

264 


0)14 3 = ALOG)DUM 2 )*AL 05 )DU«)/NP 1 

DC 13 4 

774 


IF) DU'* 3 .LT.- 8 D. )C 0 TO 1 

JD 135 

2 B 4 


1 ? P ] 1 , 3 ) = 0 A ) 3 )- )Pt ) 3 )-P()l )* 0 UM* 4 ) 1 ./MP 1 ) 

DOl 37 

294 


1 P) 1 , 3 MIP)=P) 1 , 3 ) 

00137 

304 

c 

0 RESSU 9 E PROFUt OF OTHER ABiORBING GASES 

0314 1 

314 


00 3 J- 1,32 

00144 

32 * 


JPIR 933 -J 

00145 

334 


psUh:'i . 

03146 

344 


IF)L 1 .E 0 . 1)00 TO 4 

uOliO 

354 


00 5 Lr 2 ,Ll 

0315 3 

364 


P(L, J):P) I, 3 )*CPM., 3 ) 

00154 

374 


»IL, 3 “IR)=P)L, 3 ) 

00155 

364 


4 bsU 4 =pSU 4 »p )L ,3 ) 

GDI 5 7 

394 


4 PSUMIPSU 4 ->P ) 1 , 3 ) 

60157 

40 * 

c 

RRES 5 U 9 F PROFU.F OF YON ABSORBING GASES 

..'fl 6 j 

4 1 4 


) F) L 2 .E 3 .L 1 )G 0 TO 3 

00162 

424 


0 U 4 :P 5 -PSUH 

6016 3 

434 


IF) OJH.lT . 0 . )CUM=U. 

0016 5 

44 * 


00 7 L=L 3 ,L 2 

o 0170 

454 


p)L, 3 )= 0 U«*CP )L, 3 ) 

(.3171 

464 


7 9 )L, 3 «IR)=PIL, 3 ) 

j 017 i 
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3 IF) P 5 UH. 0 T ,PS )NTE 5 T :1 

00173 

464 
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TE 4 PERATURE 30 UN 0 ARY LAYER, WALL TEMPERATURE TWB 
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cn* 


DC 6 3 = 1 , 3 B 

0077 3 
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0 U 4 = ) T) 32 1 -T) 3 b) )/)PA )I) 4 t l.-YI 3 Bl ) ) 

00205 

534 


DU 4 C:Y) 3 )/Y( 3 H) 
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544 


T) 3 ):TWn + )T) 3 p)-TWB)»OJMO-‘)DUH 4 Y) 3 B)*TWB-T)J 3 )) 4 DUMO*)DUMO-l. 

DD 2 D 7 

554 


8 T{ J 4 I 9 )=T ) 3 ) 

0020 7 

5 b* 
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?ADIA‘'CF RROFILF FOP I W A YELENGTHS , PL ANC K FUNCTION 

0021 1 

5 74 


9 00 10 3 = 1,32 

0021 4 

564 


JMI R = 33-3 

007 1 5 

594 


00 13 1 = 1 A,I 3 

60770 
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0 U 4 = 4 ?) I) /T ) J ) 

uD??l 

614 


IF)AB 5 )OU«).GT. 60 .)OUH-;PO. 

UC 223 

674 


4 ) 1 , j)=pi)i )/ )EXPinuM)-i. ) 

l.0274 
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1 ( 4 ( 1 , 34 IR)rR) 1 , 3 ) 
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r: 3 2 •' J 

654 


END 


Lnjoiz 

IjODHIZ 

ooaaiz 

ononiz 

Lnoniz 

tjDiniz 

CI700IZ 

onuoiz 

DDUOIZ 
GOOOlZ 
Lnaoiz 
uaao3t 
000037 
□ OOOMZ 
00D0S6 
LDOObl 
000063 
000066 
000075 
000107 
000113 
000. 2Z 

ODomz 

000145 

000147 

000156 

00U17D 

DDU173 

00JZD3 

000 ’03 

OCOZIO 

0C'JZZ4 

000731 

OOUZ32 

000734 

000744 

000246 

000747 

000254 

0CQ754 

600764 

000767 

000776 

000310 

onu^iD 

0OD31Z 

003316 

000316 

000327 

000344 

L0u356 

00036 1 

000366 

000371 

600406 

000406 

000412 

000420 

000423 

C0C433 

OD0435 

000443 

DD0452 

006462 

LDU''23 
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APPENDIX F 

TRANS SUBROUTINE LISTING 


!.^ni 


i-Olo ( 


L.ol'd 


3010 , 

M * 

roi.' J 


30103 


301" 3 



0# 

^.OlOb 


0 JlDo 

n-t- 

0010 7 

ii<f 

1.010 7 

iSfl' 

001D7 

13« 

00107 

14« 

LOi 1 

IS# 

0012O 

1 ft# 

00123 

17# 

00126 

1ft# 

0C1?6 

19# 

0312 1 

?7# 

00172 

21# 

001 39 

22# 

COl 3x 

2’# 

JOI 35 

24# 

0013 7 

25* 

3019 1 

2(3# 

03 l‘<2 

27# 

Dr 199 

# 

i;oi95 

?^# 

JC19 0 

’I r# 

0019 7 

n# 

tiOlB.. 

32* 

r 015b 

37 # 

,0157 

?4# 

0015 7 

35* 

LOlbu 

3ft* 

u 0 1 5 t 

37 * 

0015m 

3P# 

03165 

39 * 

00167 

4a* 

V.C167 

4 1# 

00173 

42* 

00176 

4 3* 

,.317 7 

44* 

00202 

4 S# 

00203 

4ft# 

05236 

47* 

0L21L 

4ft* 

UD21 3 

49# 

002 19 

5''* 

0021 X 

51# 

00220 

52* 

002 2 3 

S3* 

03226 

54# 

0022 7 

55* 

00270 

5ft# 

00231 

57# 

0C232 

ftft* 

002 3 5 

59# 

1.0236 

ft'^# 

E:0'2 3 7 

61# 

002 9 0 

ft?# 

0029 1 

ft3# 

1.D292 

^4# 

00392 

65# 

0329 6 

66* 

00 251 

67* 

E.026X 

bft# 

COl'Xb 

69* 

6025 b 

7*1# 

0025 7 

71# 


^UHRDimvr 

11 , j:, ,psuM,jB,i?i(9),9?(Ei>,x(Hn, 

X P«(4),ri4uUR),Ci;L(31,T(Ul|,'iximi,><X,Nti,!>5,»ATH,L?, V(?n,L3, 

X TXTfc'M P) ,B3AM|B,3,R) ,K 0 (A, 3, 211 ,kK( 9,3,ul ) ,riTtB (8) , 

X lfA8(8,3,B|,SB(6,3,?ll,sR(R,-',m),BU8,3,m),ETA(B,3,<(l), 

X U(R,’,All,TAU(P,A,Al),NPl,rA,IB,MT!:ST,TX 
PEAL KX ,<TAB,MB ,NP1 
I'jrrsrp ppop 

nr ON BBETA"( 7) , A1 ( 7) , AZ|7T,A3( 7| ,AAI 71 

5ATA '’BEnB/D.,?.,B.,ia.,?a.,53.,103./,Al/-.68,-.37, .AB,3.1P,9.13, 

X ?9.1,lPb./,4B/-.3E,.(,B,3.1i,P.25,22.n,B7.1,239./,A3/.56,.71,.9P, 

X 1. lb, 1.67, 2. 76, A. 35/, AA/, 07, .A3, .79, 1.33, 2. 26, A. ?6, 7. C9/ 

C BA^O “ortl. PAPA^ETEA PPOXILE "7 IABUlAR I A TEPPOLA TJ ON 

no I 

30 1 in, El 
DO 2 J=1 , J? 

JPIBPJ3-J 

C interpolate FOS KL and 9D, OENOPAALIZC KD to XX, EOllA.l) 

DC 3 N=2,NK 

If ( T ( Jl ,BT .TKl AB 1 M.ANn.N.NE.NrtlGC TO 3 
•<Cn,L,J):KTAn(T,L,N-ll*IKTAB(l,E,N)-KTAri(I,L,N-in 
X»(T(J)-TkTAPI N- 1 M/ITXTAL'( N)-tkTAA( N-in 
If(T(J!..T.TXTAanMKUl,L,J)rKTAR(T,L,l) 
IF(T(J).3 T.TkTAR(NH)|K3(I,L,J|:KIA8(I,L,NX1 
30 TD A 
: OCNTINIIF 

A <X(I,L,J):aC(I,L,J)*273,ap(L,J|/T(J> 

2 <k(IiE,JAIP|:kkII,L,JI 
c interpolate fop tt'.. 

DO 1 J ' , J2 
DC 5 "A 

l'( Tl J) .bT .T jTAR ( N) .AND.N.NE.NbI 60 To 5 

a;ii,l,J)patae< (i,l,n-ii*(ptao(i,l,n1'"tabii,l,n-i)i 

XA(I< Jl-TbTABf N-I ) )/(I0TAP( M-tbTaBI N-Hl 
IK T( J) .uT.TATAP I 1) )3U r,L, J|:BTAH( r ,L ,1) 
IFITIJI.CT.TKA9('.K||j3(I,L,JUPTAB(l,L,NX) 

3 0 TO 1 
b CONTlNUr 
1 oontinuf 

C -INE PPDARENINC AND DEN0ARAL17E 30 TO BP E0(1A.2I 

DC 6 J-1 , J2 
3-J 

DO 6 L:I,L1 
DUMRD. 

30 7 LLP1,L2 

7 DUM:DUM + ALPH|L,LU»P(LL,JI 
DO 6 T=TA,IB 
3B( I ,L , J ) rBOl I ,L , J)*DJM 
6 3b( I ,L, JHIPl-BBI I ,L , J) 

C 3E E3(l?> 

DO S I=IA,ia 

DO B L::l,Ll 

3E( I , L, n TAB 1 1 ,L , 1) 

S1:D. 

52 = 0. 

SAVlPXiiU ,L, 1 )*PB(1 ,L,n 
DC » J:7,J1 

5 AI/2:XK I I ,L , J I7PB ( I ,L , J I 

51:51 + 1XMI,L,J-1)+KX(I,L,JI)P0X(J) 

52:S2*(SAVl4£AA2)*DX( J1 
3AV1:5AV2 
BE< I ,L,J>TS2/S1 
A C0N71NUF 

C OPTICAL DRPTH U EOdll, AND ETa -OSinddA) 

DO 9 I =I A,IB 
DO 9 L:I ,L1 
E TA ( I ,L , I ):C. 

Jd ,L, d:'. 

SAVdKK ( I ,L, 1 ) ORIGINAL PAGE IS 

OF POOR QUALITY 


0013126 

003126 

003126 

000126 

000126 

003126 

000126 

OOUUb 

00O126 

U00126 

003126 
00U126 
000126 
000126 
003126 
00013S 
000166 
003 166 
000166 
000175 
000175 
000212 
000212 
00022A 
ODD 2x0 
ono’50 
00u?b5 
DP02S5 
00Q261 
000261 
000^7 1 
000311 

uOOTd 
000326 
00j''26 
OOU3AC 
COD 3A6 
C0035A 
00Q376 
003 776 
000376 
one 376 
DOO 376 
OOOAUl 
LODUIO 
000912 
OnuA12 
0P3A2C 
000920 
0009... 
000922 
00UA56 
000956 
OOOABt. 
000957 
000960 
C0Q9E,1 
00 3 9 b 6 
000966 
0P097D 
000975 
OnORCl 
000R03 
00u523 
LDur23 
00052 3 
003532 
003537 
000590 
00059 1 
000592 
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00260 

72* 


DO 9 J=2,J1 

000551 

00263 

7 3* 


385:88(1, L,J»/BriItO,J) 

D0D551 

0D,’6‘) 

74* 


I F( BBO.nr . 100. |PBC:100. 

000554 

Q0266 

7 5* 


SAV2:KK( 1 ,0 , J»*PBi:##5T« (I ,L,J-1 ) 

000562 

00267 

76* 


Sl:Sl*(SAVl*S6V2t*DX(J)#.5 

OOu572 

00270 

77* 


SAV1:SAV2 

000577 

00271 

78* 


J(I,L,J):S1/BE(I,L,JI 

000601 

00272 

79* 


DO ID S:2,7 

QI306U6 

00276 

80* 


IF(BBE.GT.3SCTAP(tj|.AND.»i.NE.7IG0 TO ID 

000606 

00277 

81* 


9ATI 0:( 3BE-BBETAb(*J-l ) )/( BBET Aflni -8faEIAB( N-n 1 

000623 

003QQ 

82* 


A5:l.+.ia5#B3E 

0DD631 

00301 

83* 


IF(U(I,0,J>.GT.A5IETAn,0,J):(U(I,L,J>*AllN-', •♦(AK'n-AKN-lIl 

000635 

00301 

84* 


X*BATIOI/(U(:,L,JI*(A2(Nl-A2(N-lH*RAiro*A2(»J-lll 

000635 

00307 

85* 


IF(UII,L,J).LE.A5)ETAn,0,J):(Ua,L,J)#IA3IM-ll*|A3(N)-A3<N-in 

000656 

00303 

86* 


X*RATIOn/(A5 + U(I,0,J)*((A4(NI-A4IM-n)*RATIO*A4(»J-lin 

00D656 

00303 

87* 


GO TO 11 

DCD677 

00306 

88# 

10 

:o>niNUE 

000726 

U03I0 

89* 

1 1 

C04TINUF 

0Cj72S 

0031 1 

90* 

9 

OO^TI RUE 

00U726 

0031 1 

91# 

0 

T AO FOR L 1 GASES, EO (9 ) , (13 1 

0DD726 

0031 6 

92* 


DO 12 l:l#,13 

000726 

00320 

93* 


DO 13 J:l , J1 

000*26 

00323 

94* 


TAUCI ,4,Jl:l. 

0DD725 

00324 

9 5* 


DO 14 0:1 ,01 

000731 

00 32 7 

96# 


T A0( I ,0, JUEXP ( -U (I ,0 , J)#BE(I ,0, J) /SC)HT( 1 . *1 .570B*U(1 ,0 , J1 n 

000731 

0033b 

97* 

14 

TAU(I,4,JI=TAU(I,4,J)*TAU(I,L,J) 

00075 1 

00332 

98# 

13 

OORTlNUr 

000766 

U0334 

99* 

12 

OTAULdCTAUn ,4,J1 

000766 

03334 

133* 

C 

3AS RADIARCE CRL E0(8) 

000766 

00336 

101* 


DO 15 I:IA , IB 

001015 

0034 1 

102* 


CRL(Il:l3. 

001017 

00342 

133* 


DO 15 J:2,J1 

001"22 

00345 

104* 


DTAU( I ,J|:TAO (I ,4, J-1 )-7A0II ,4 , J) 

001022 

03346 

135* 

15 

ORK n :C50 (I 1 • ( R ( 1 , J 1 *R ( I , J- 1 M #DT AU ( I , J 1 # . 5 / ( 1 . -0 T A 00 ( 1 1 1 

CD1024 

U03S1 

106# 


RET0R4 

jDl'’45 

00352 

107# 


END 

001 122 
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APPENDIX G 
WAVES PROGRAM 
Program Symbols 


The following FORTRAN symbol list is limited to symbols not in the 
PROFIL symbol list, APPENDIX D, and not in the NAMELIST data in the 
section WAVES PROGRAM USE. 


FORTRAN SYMBOL 
AMAX 

AMIN 

KEY 

KEY4 (5) 

K5 

NPB 

PASAV (8,10) 

PASAVE (4) 

ZETA (11) 


DESCRIPTION 

Maximum value of dependent profile parameter to 
determine ZETA 

Minimum value of dependent profile parameter to 
determine ZETA 

Selector for variable profile parameters^ KEY = k4+ 
PROF-1 

Index parameter (data statement) 

Index parameter to select dependent variable profile 
parameter 

Index parameter NPA-1 

Dependent variable parameter after correction by 
ZETA 

Saves initial value of independent variable parameter 
Graph -plotting constant (APPENDIX C) 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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APPENDIX G 
WAVES PROGRAM 


00101 

1* 

0013 1 

2* 

00101 

3* 

00101 

4* 

00101 

5* 

00101 

6» 

00103 

7> 

0E1D4 

BA 

03136 

9* 

UOlDb 

lOA 

00106 

11* 

UU10 6 

12* 

COlDb 

13* 

00137 

14* 

0010 7 

15* 

0010 7 

16* 

03113 

17* 

LOllO 

IB* 

03113 

19* 

00121 

20* 

U0122 

21* 

03126 

ZZ* 

00127 

23* 

001 32 

24* 

0014 3 

Z5* 

00157 

26* 

00173 

27* 

0017b 

28* 

00177 

29* 

03232 

33* 

00206 

31* 

00210 

32* 

00213 

33* 

00220 

34* 

00222 

35* 

00226 

36* 

00230 

3 7* 

U0Z41 

38* 

00256 

39* 

00255 

40* 

U0271 

41* 

00273 

42* 

00275 

43* 

00276 

44* 


COMMON P!!0F,TWb, LI, Jl, JZiJ’iCl .PSUM,JiSRll9 ) ,■<:(« I ,M41 I, 

X P)^(^),CTtUUI9),CRL(B^,T(4n,nx(l^ll,'|K,^l3,PS,P^IH,L^, Y<Z1),L3, 

X CTAU(8,4n,Ptb,41l,CP(6,Mll,Pt6,4l)>''LOH(J,B), 

X TXJAP( 6) ,KTAB (B , 3,P) ,KD(P, 3, 2U ,XM 9 , 3,41) .TPTftB (8) , 

X RTAB(P,3,e),6Cl)S,3,21),8B(8,3,4n,PE<8,3,4n,ETM8,3,41), 

X U(B,3,41l,TAU(8,4,41),4Pl,lA,IB,NTrST 
SEAL KK,KTAB,HB,NP1 
INTEBES POOF 

OIHEN8I04 DPAPAl4),PA9t8l,DPA(4),nELI4),RA(8),[JPV(3),PAlll0)i 
X B0TLl4(4l,UPUlP14|,HRL(B),MTAUL|P),GI9,b),GV)9,b),GAS|(j),UAVE(4), 

X Z( 21 ) ,0b ) 8 ) ,PASAV( 8 , 10 1 , P ASA VE ( 4 1 , M A4£ ( M I , 2 LT A ( 1 1 ) , KEY 1 I 5 ) , 

X KEYZI5),KEY3(5l,KEY4(5l,PUtb),NP(6),PC(6| 

seal MrAUL,MSL,M ,NP,HAX 

MAKE LIS I /N AMI /IRE AD, POOF, PATH, 11 ,WAVF, 

X TW,TC,4T,TKB,JP,UPL1H,ROTL1H,L1,L2,PS,PC,1P,Ph,».K, 

XIK7AB,9P,TB7AB,H4,NPA,PA1,J1,Z,ITEH 
DA7A nEL/4*Ol/,NAHE/ZHN7,4HTC,K,6HPC,A7M,4M7W,K/,JI/7l/,?/0.,.ni, 
X.DZ,.r5,.l,,Z,.7,.4,.5,.7,l.,I3Pl*/,XEYl/Z,l,3,D,?/,KFYZ/3,3,4, 

X 0,3/,KEY3/l,2,l,n,l/,KEY4/Z,I,4,D,?/ 

600 IREAD=n 

READ! S ,MA«1 ,END:999 ) 

IF( IRrAn.FO.DIGO 70 113 
RCA018,BD) IGASai ,L = 1,6) 

READI S,Pl ) C 1 ALPHIL,LL ) ,LL = 1 ,8 1 ,L= 1,L 1 I 

RE40(5,B11l(IKTABn,L,N),N:l,P),IRl,Il),L:l,Ll) 

REAO(5,BllIII6TABI7,L,N),NPl,P),I=l,Ill,UPl,Ll) 

00 bb L = 1 ,L1 
ALPH(L,L)=K 
DO bb l:l ,I 1 
00 bB N:1,NK 

68 IF(KTAB(I,L,N).L7. . 1 E -3 6 ) K 7 AB <1 ,L , N ) - . IE -3 b 
DO bb N = 1 ,MB 

86 1 F( BT AB( 1 ,L ,M ,L7 . . 1 E -3 b ) B 7A3 U , L , N 1 = . lE-06 
IID HRI7E16,9D) 

i4RI7E(b,'IA‘n) 

WRI7E(8,36> <GASIL),L:l,b) 

KBI7E(6,b7)((ALPH(L,LL),LL91,F|,L:l,LII 
ilRnE(b,BB)((IKTA9tl,L,N»,N: 1,81, 1=1, III, L = 1,L11 
WRITE(6,B9»((|B7AB(I,U,N|,N:1,8),I=1,I1I,L:1,L1) 

C INI7IAL CO'JDllIONS 

IFIK4.E0. 1 IPROFrl 
I F( X4 .E0.4 )PROF=Z 
<EY=K4«PR0F-1 
XlRKEYl UEYl 


CDOPOD 
LOoroD 
3D3C3D 
DDurUD 
GCjPOD 
uOLOLO 
LPO"Lil 
LCunoi 
uDunol 
oDoroi 
ooLnoi 
LDjnD 1 
Drug’Ll 
DDUPOl 
tnonoi 
LOUPCl 
LCLPLZ 
COJBOZ 
LDuDuZ 
LDjrOZ 
EOGDEZ 
DD0DU7 
cnoDi 1 
CDJP43 

CC3113 

LDU144 

000144 

LCalbO 

000160 

003160 

OOOlbl 

000161 

00D7U0 

ono?D4 

cnoRio 

0DD756 
000302 
L0J326 
000326 
OOU353 
000360 
000365 
00037 1 
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00277 

45* 

<2:iteV2(<eVl 

00300 

46* 

<3SK6VSMevi 

00301 

4 76 

<5:K6V4|HeVI 

00302 

466 

9TES7:0 

00393 

466 

l«Sl 

00304 

506 

I6SI1 

00305 

516 

J3SJ161 

00336 

62* 

J2s<Jl6||/2 

00307 

536 

b3Sbl6l 

03313 

546 

eispuni/oeiii 

00311 

656 

66I1ISN7 

00312 

566 

P*l 2 t:Te 

00313 

576 

6«l3i:6C(ll 

00314 

SS6 

6A(4I:tw 

09316 

896 

00 64 Kslt4 

00320 

696 

9P*6*IHl:0. 

U0321 

616 

DPAMISO* 

00322 

626 

PA66VEIK>:P6|H) 

00323 

636 

90 64 |:it6 

00326 

646 

BIKtllsO. 

00327 

656 

01864, l|:n. 

00330 

666 

Sl».l|:9. 

00331 

676 

»*S*VIK,l6«i:o. 

00332 

666 

6«SAVIK64,l64|:0. 

00333 

696 

00 64 4:i,41 

03336 

736 

7AUII62,h,M1:C. 

Q0S3T 

7|6 

*11,8150. 

00343 

726 

CPtl,P|59. 

00341 

736 

64 C097I4UE 

09346 

746 

90 59 |:i,ll 

00350 

756 

411115 .ll909EoU 6UAVEI1I663 

00351 

766 

50 92III5l.43BB6WAVEm 

00353 

776 

90 61 JS1.J2 

110356 

766 

J8IRSJ3-J 

09357 

796 

«|j|5P«T862U)/2. 

0036b 

606 

V1J152IJI 

U0361 

616 

51 8IJHIRI5PATH-XIJI 

00363 

626 

90 59 J52.J1 

00366 

836 

59 9XIJ|5X(Jt-)(|J.n 

00366 

646 

C PBESSOSr PROFILE INtllAt PROFILES 

00370 

BS6 

4P1S4PIII 

1)0571 

866 

90 52 J51.J2 

00376 

876 

90 S3 LS1.L2 

00377 

686 

3U81:A6SIPCIL >*puILM 

00400 

696 

P|L,JI5PCILI 

00401 

906 

IFIDUH1.CC,0,.0R.J.E0.J2IS0 TO S3 

00403 

916 

OUHSl.'Vt Jt 

>.0404 

926 

OU825AL06IOUPn*ALOSlDliPI/NPILI 

00406 

9 36 

IFI0UM2.iT. -69.160 to 53 

00407 

946 

66 PlL,JlsBClL»-lPClLI-P«llLn60UM66| 

004 lU 

966 

53 COMTIVUF 

00412 

966 

IFILI.EO.IIGO Tfl 55 

09412 

9 76 

C PRESSURE PROFILE RATIO FACTORS CP 

00414 

986 

90 54 L52,L1 

US417 

996 

54 CPIL,U»5P(L,JI/P11,J> 

00421 

1006 

65 tFlL2.E0.Lll6O TO 52 

U0423 

’ 16 

DUM50. 


003 174 
000776 
000400 
000402 
OOOttOS 
003400 
000407 
0004 It 
000413 
000416 
000420 
000422 
000424 
000426 
0004 65 
000465 
000465 
003456 
000464 
000464 
000464 
000466 
000466 
000467 
000472 
000472 
000472 
000473 
OOOS17 
C0U4I7 
003617 
DD0S23 
00U632 
0PP»J2 
000636 
000661 
000643 
ODOS63 
000653 
000663 
0006S6 
000666 
000606 
000606 
000610 
000612 
000517 
000621 
000633 
000636 
000666 
000656 
000666 
000660 
000666 
003672 
000674 


ORIGINAL PAGS IS 
OP POOR QUALmr 
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00U24 

C0427 

00431 

00434 

00436 

00440 

00441 

00442 
00446 
00446 
00450 
C0464 
D0466 
00466 
00470 

00473 

00474 

00475 

00533 

00501 
00532 

00502 

00503 

00534 
00505 

00510 
bOSIO 

005 11 

00512 

00513 

00514 

00515 

00515 

00516 
00530 
00523 

00521 

00522 
03535 

00536 

00537 

00530 

00531 

00532 
005 S3 

00534 

00535 
00636 
00537 
30540 

00541 

00542 

00544 

00545 

00545 

00546 
0C547 


102* 


30 56 I.:i7,t.2 

00U7U0 

133* 


56 9UHs0UM*Pa,J> 

000710 

104* 


30 57 LSUSiU 

000717 

105* 


57 ePHiJ»s*a*JI/0UM 

C0Q7I7 

136* 


52 eOVTtOUE 

000726 

107* 


e«lL TE*P 

OOOT26 

lOQ* 


CtLl TP64S 

000730 

109* 


30 65 tsi.n 

000732 

no* 


4ALIIISC<UXI 

000740 

111* 


63 *muimsC76UUTl 

000741 

112* 


*AITei6»3SlHC*$UI<L:U6nmjn7iJi,tPn,JI,tsi,6I.J:t»jn 
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DO 22 KSKl.H2.H3 
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test fOR RBRBMEtER OUtSlOE ItHIt 
DO 23 KS1.4 
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35X FORHATIIHO.IPMTEHPERATuRE profile, SX.3BHPARTIAL PRESSURE PROFILES 
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EXAMPLE or PROFILE PROCUIAU DATA CARD! AND OUTPUT IN MODE - I 
Dtti Cirdi C^NniMd Prom Oirdi of APPENDDt % Wllli DtM Pbr 
Pour Wivo Numbori <tl -4), ProOle moP • I 


SNAH1 M00E*3«n*4 $ 
W PN*.»*.<4»4*0. 


$ 


Printout of Uipttt Dnin 


«n*Ht 


HOOC 

s 

♦ 1 




n»of 

s 

*2 




i»e*o 

s. 

*3 




l**TH 

s 

,toooaoooc*3s 




n 

s 

♦* 




HAVE 

s 

.S$939ec9e*34i 

•310909091 *34 • 

•349399391*34. 

.410000001*04. 



•50003000e«3B| 

.000000001*00. 

.000990991*90. 

.000000901*00 

HH 

s 

•3760»72«E*92» 

•3772B00tl>02. 

•37I640741-02. 

. 3 B 71 * 27 lt*n 2 . 



»37332n37C-92, 

.000300091*33. 

.093399991*90. 

.999033991*09 

HTtUL 

5 

.99S3I32BE*00» 

.822343711*00. 

.802994tSe*00, 

• l6S6i»SSt*00. 



.BS4I33SBE*33. 

.090900D01*00. 

•000000901*00. 

.000009001*00 

&!?«« 

s 

.IOOQ9COOe*91 




e»tn 

s 

•BQ0aaB09E*01 




ru 

z 

*17393909£*9« 




tc 

s 

.200090001*94 




NT 

s 

•300090001*00 




lut 

z 

•909330091*99 




UB 

s 

*1 




unuM 

z 

. 100090091*01 • 

.100000001*04. 

.103990991*01. 

•300000001*04 

BOTLtH 

z 

.2a9390C9c-9l> 

.690909091*33. 

.29933999101. 

,609039301*03 

U 

s 

*1 




L2 

z 

*2 




PS 

z 

•199930091*91 




PC 

z 

.3S9D30D0E*30. 

.630000001*00. 

.009990391*90, 

.000099991*00, 



.00003COQC*90f 

•000000001*00 



NP 

z 

.S0399(<t3C*93. 

.390309091*93. 

•039399391*90. 

•993039301*00. 



.900090091*901 

•000000091*00 



PW 

5 

.300090001 430t 

.700000001*00. 

.000000901*90. 

.000000001*00, 



.90O03nO0E*90» 

.030909091*93 



NN 

z 

*» 




nun 

z 

.B0a03Q00C*33f 

•100000001*04. 

.169390991*04, 

.209000001*04, 



.2S9339t'9E*3«, 

.330399991*94, 

.09)333391*30. 

•393099391*09 

NB 

s 

*B 




IBT4P 

z 

•BC0090C0E*3S» 

.100000001*04. 

•164000301*94. 

• 200U00e0E*O<i, 



•2S0990001*34i 

.390309091*94. 

.a393)933E«33, 

.393033331*09 

J1 

z 

<21 




2 

m 

•39939Db3E*39« 

.looaoooaeoit 

.2333)g30C>01, 

•600000001 *01. 



•t00990D0E*30t 

.200000001*90. 

.309390901*90. 

.400000001*00. 



•300090001*90. 

.700000001*00. 

•100390391*01. 

.109000001*01, 



.ia93900ac*31» 

.100900001*91. 

.199339301*31, 

•109000901*01. 



.100030091*01. 

.100900001*01. 

.103390391*01. 

•100000001*01. 



.193330001*31 




ITES 

z 

*10 




It 

• 

.900090001*00 




lENO 
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M 09 C : 

Pfor s 

Item s 

PtfH t 

II i 

tt«ve I 

H«k t 

Hflui I 

{••v X 

CKth t 

1M 1 

1C X 

HI S 

TMI X 

Jl X 

UHin X 

■OTkIH X 

LI X 

L 2 : 

PI 

PC s 

K* X 

PW X 

NH X 

tHT*P I 

hi X 

riT«n • 

Jl X 

i X 


ITEI 

TX 


*2 


♦ 3 

•IOO 930 C 3 r *31 

• I 

.110030031 *31 • 

•Se003CD3E*3X. 

•I7»3*T2IE*3S. 

.}73S23in>32. 

*IM}|S2ie*30i 

.•»«l33IIE*33f 

*I0093CC3C*OI 

•133330(132*31 

•ITa330032*3« 

• 2D0S30C82*3i> 
*3333)303203 
•30003DCaE*30 


*1 

*l3333903e*)|« 

•2C0030032>3t. 

»l 

♦ 2 

•103330032*31 

•)K3333C3E*38f 

•)33333r32^33» 

•5CaC3?C32*30f 

•00003r032*38« 

.303}3r3e*)3. 

•3CQ030C3C*90i 

•I 

•103033032*33. 
• 2!033r<C32*3». 

»k 

•100030032*33. 

•22D333C32*?». 

*21 

•3cooaooor*33. 

•103030002*30. 

•233339232*33. 

*100033032*91. 

.193330102*31. 

.100030032*31 

♦10 

.333333232*33 


• )IQDC0O0C*0<i. 
.000000002*00. 
.ir72l99f2>)2. 
.900000002*09. 
•l223l«7ie*D0. 
•000000002*03. 


•139333032*31. 
• <>00000002 *03. 


•^«CJ0Cn02*33. 

•990303031*33 

•200000002*00. 

•000000902*03 

.&«03S9932*33. 

• ''00000902*00 

•I00DDD002*0». 

.^000009or*Qit. 

.t0ODODO02*Oti. 

•100000002*01. 

•lOoceoFOE'Oi. 

•?00000002*00. 

.733333332*3). 

•100000002*01. 

• 100000002*01 . 


• I«a39039i*3«. 
•003*30332*00. 
•)7S2tl7«C*)2. 
.000*00392*00. 
.M29mS2*39. 
•003330332*30. 


.133133)02*91. 

•203933331.31, 


•033139)32*00. 
•003*00392*30. 
• 0339)3)32*3!*. 


• )33933332*3I». 
.003933312*00. 

• ISU90030t*0«i 
•093930332*00, 

• 20u90030L.31 • 
•100330302*00. 
.13)333332*31. 
•103900302*31. 

• 100000302*31 . 


12h3 


ORIGmL PAGB IS 
OF POOR aUAUlY 


.4IOO0393f*rit, 

.OI)aODo002*ou 

.J*71f37Sf.!i*, 

•000000092*00 

•n}i«»)i2*oo, 

•033033332*93 


• 3ooooo'*oc*nit 

.603003902*03 


.903003092.90, 
.2103003002 *00, 
.33303393E.93, 


•233033*32*96. 

.009030902*00 

.200093932*2 4, 
.30300#*32*9C 

#S0000i)ODC-9J, 

»4C0000O02»*0, 

•133033132*91, 

.100000902*01, 

.100000002**1, 


CMKWUtton and Fnifltti 


P'O^IU 

ff pp 


turn 

HIU 


PtK'S 


t»»OKCVt 

te-*? 

p«p»s 

n>*p 

PBCSS 

r*»e*(tNr 


Vt 

K,it 


t4«H 




>f>(.33S 


.H-IPO 

tP33>D3 

•IkBSO 




Jpt. 

,33tt,3 

.'J 

.B3"3a 

♦33 

)T£RI?I04I 

3 94$ «e$tCU6bs 

.394-31 

$tep* .009 

p»#*»eTE»i 

n 

•«3SD3 

3U3P.C3 

.S»"CiO 

1T33.9C »£ttPUII.$: 

lilVthimiEPU/CMlt 

•♦Ul-"l -.749-0? - 

llwO. llCOt 

.105-31 

3400. 

-.692-02 ,S49«r.i 

4|00. 5500. 

lNC-lt*1tHTS 


.33RC3 

).lb 

•.0»»2<* 

13.42 

ncRiriovs 

1 945 «E$|p-3*U 

.150.32 

5?ePS ,129-01 


PI 

•R33L3 

2331,t3 

.33IT6 

1313.42 P:5IQUII.5i 

lIVENliHBCPU/CHl: 

•.629-oS ,44?-04 • 

35wfl« 3tr,o. 

.631-33 

34D0. 

.620-03 ,117-02 

4100, 35or. 

iNc»r«vTs 

3P* 

• 33.U3 

•1.51 

-.33353 

1.49 

MERIU 0N= 

2 445 n;5IDUII.s 

•626-03 

5UPS .311-62 

P*»»BETt»» 

PI 

•I33C3 

?03t.?» 

.331*2 

1713.91 PEilDUIkks 

wivehUHecP(t/eHti 

-.4)0-95 -.372-03 - 

3509. 3109. 

.226-33 

3400. 

.445*03 594-04 

41b0, 3^00, 

lNB«rHE«ITI 

3PI 

.D3QC3 

-.Pb 

(DOPbl 

.13 

ITCRIItO^i 

3 9M$ RE5IDUII.S 

.62 5-03 

SIFPS ,ll3-b3 

»4»*M£TE»S 

PI 

.S3303 

2331.33 

.33125 

1714*03 PE510UII.5: 

Hive hUPPC*>ll/rPli 

♦2BP-"4 375-33 ' 

35)3. 3103. 

.294-33 

34(i0. 

,463-03 -.709-06 

4)b0. 35(.r, 

iNC^EdrMS 

3 PI 

.33313 

.30 

.CbCbO 

.00 

ITEPIllONS 

4 94$ 9C5IUUIIS 

.623-93 

SIEPS . 929-05 

p»»*METe»s 

Pi 

.5331(3 

2331.23 

.33126 

1714.54 PtSIOOItSs 

4IV(.NUHr<LPt)/(M|S 

,2»7-54 -,375-33 • 

35uO. 31)9, 

.233-33 

3400. 

.«5!-1)3 -.521-05 

4100, 3500, 

nu»tH« caooR 
PCRCt'lT 

RMS TOTtl 

.31 

.3 

.3 

.3 

.3 

.0 

-t.lt 

?.i( 

*.T 

t.) 

*.i 

1*5 

.1.5 

-.t 

.2 

-.5 

T.l 

t.l 

« 

3.25 

•4.5 

1.9 

•1.9 

.2 

6*6 

6VEHJ46E4I1/C4I 9£UnUH.PeRCE4T 
3530. 7.0« 

3130. 1.2 

3400. 2.0 

41)9* 3.6 

3500. 4.4 

0. .0 



tUVENUMREBn/C^I 

f»*mSHTTT*»»CC 

s JE33. 

: .5)1 

3133. 

.«22 

3433. 

.P02 

413). 

.153 

3)33. 3. 

.0)0 .009 

3. 3. 

.000 .000 




TtMPEB*TU»C Pftorlie BBlSSUUt PBOFIltS ClfH St&vATH 

*,cp rcMPi)i Hio s? 


.33 

1714,34 

• 3437 

.6593 

• 0009 

.OOuO 

.”300 

.0000 

,50 

1719.75 

.3405 

.5595 

.0900 

.0003 

.3300 

.0300 

l.OC 

1725.41 

.3433 

.5597 

.0003 

.93)3 

.3330 

.0303 

2.50 

1742.94 

.3595 

.5602 

• 0000 

.oobo 

.3000 

.0000 

5.00 

1760.63 

.3389 

.6611 

.3333 

.0033 

• 3300 

• 0303 

tb.oo 

1617.43 

.3373 

.5627 

.0009 

.0033 

.3360 

• 0000 

15.33 

1653,51 

.3359 

.6641 

>3333 

.0033 

.3300 

• 0309 

2b. QO 

1197.64 

.3347 

• 5653 

.0009 

.0003 

.3300 

• 0900 

25,33 

1929.43 

.3336 

.6654 

.0009 

.0000 

• 0330 

.9000 

35.00 

1975.36 

.332) 

.5679 

.0000 

.0003 

.3300 

.0900 

5b.00 

2901.23 

.3313 

.6667 

• 0003 

.0000 

• 0300 

.0900 

65.00 

1975.36 

.3321 

.6679 

.3333 

• 0033 

•9300 

.3303 

75.00 

1929.43 

.3335 

• 6664 

.3000 

.0003 

.3000 

.0900 

6J.33 

1697,64 

.3347 

.6653 

.3333 

.0333 

• 0900 

.0000 

55.00 

3660.5) 

.3359 

*6641 

.0000 

.0090 

.0300 

.0000 

93.30 

1617.43 

.3373 

.6627 

.0000 

.000’* 

.0900 

.0000 

95.00 

1766.60 

.3369 

.6611 

.0333 

• 0333 

.9390 

.3303 

97.50 

1742.94 

.3391 

.6602 

.0000 

•0033 

.9000 

.0300 

99.00 

1725.41 

.3403 

.6597 

.3333 

.0333 

.9390 

.3303 

99.10 

1719.75 

.3405 

.5495 

• 3003 

.0333 

.9300 

.0000 

133.33 

1714.34 , 

.3437 

.6693 

.0000 

.0000 

.0000 

.0000 
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»!f 0 N? 


•ii» .«]| 


”Ftv=i H; 




•?*» •Jii 
:*»« tiH 

‘Osaj 

»JJ4 ,413 

•?s* •*»» 

*J 7 ? -nj 

lyij 


•ig» 
*JH0 ,>)?to 
f»n i»,o 

KJ6 g,tf 

••M «.2« 


P*or - M 

p*Th - ♦! 

ti - ••OOOD0£i3t*OJ 

■AV^ * *A 

.JiooasaoEo*. 
?u : •“SJJ^ecsrop, 

»C i •*M539f'3£«3(, 

nr . 'JposJODaEto# 

»«s - *:!??*oooeoo 

gd j •3333J5rJE»J3 

udtlM ’ «i 

soTtm : 

11 - •*M 330 C 3 E. 3 j 

»•? ! *1 

PC j ’J?£??°®3E»31 

N» . •f53339C3E*33, 

•®°®3D8CBr*3oJ 

P« . ’?P 9 ? 30 C 3 C« 3 o: 

•133333C3E.33; 

NH ; ‘“COOSOOBE^OO* 

♦* 

•t0B330D3£O3, 
yfl . •<60B3000E*34, 

* •*I^003000E»33, 

<1 . •2SaB3S03£o». 

•I : ♦! 

«| . *4 

•60003DOOE-11. 

•333335t!3^33, 

1 . •60003qo3£«Bo! 1 

- ♦Bl 

*?S£2?®°®3*33, , 

•*SS 33003 E. 3 o* ! 

•S0333CtOE*33, . 

tl0033303E*3i. 

• j 00 B 3 n 03 E« 3 j, * 

pQ . • i O0030OQP ^01 “ *' 

♦iO 


Pr/atoutoftqiut I 


• 550D0Cf'0r*o<i 


• iSMogo?**?!' *^l5glgoggt,jj 


•ftSOOOBnuE*33. 

•^303B380E«33 

•^Joooonor^oB, 

•JJpOOOfOCfQs’ 

•^303 boD 3E*33, 

fOOOBCBPOE»08 

• 1008B080E»B«. 
•3BCOOOPOC»0»I 

*jOOOOOBOE»o». 

•iacaooosE» 3 *J 


•IOOOOOPO£*03, 

••333B3B3E*33,' 

. 6 e 000 B 0 CC »00 

•JOOBBOnar. 31 , 

»?0P8DCn0£*BB, 

•’ 00300 floc» 33 ! 

•iBoaoono£»oi. 

» looBOBoactsi , 


•Oa333333E*3C, 

•O0ai0939£»3o^ 

•03333333£«3p, 

•'*353333E*3», 

•O 0333 O 39 E.O 0 , 

•*5353033 £»oi,, 

•^^?5»333E.3e, 

• S33J.'333£*3B, 


•ao 303 aPBc»ao, 

•OBoopaaoE^no, 

• 3 a 3 aS 333 E*PB, 

•£33e3ai3to<», 

»ooofloanof* 9 B 

•f“Soa303E*n4. 

•I0000300E»(JO 

• 2 BOBOOOO£*nB, 

•63B33a33E»3a, 


• Joajjlosjl^oJ' ’'ososJisr.Bi, 
•»a 33 :i 333 E* 3 ?,’ 'JS 3 S 3 “»of * 0 ^ 

• «oanti 039 E* 0 i, 

•10BB3B33f»3l, '?££££°®®'*aj » 

■”» •^OBBBOOOr»01, 
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Profiles and Computation 


t£M>e»«TO*E PUOritE PP£SSU»£ PBOrUES E»CH S*S,*TH 

M«CM rrHPf' ^2^ N2 


fOO 

1300.30 

.3303 

.7 000 

• 0000 

.0033 

.0300 

.0000 

.S3 

1)39.93 

.3)13 

.6990 

.0000 

• Qoon 

• 0300 

.0000 

1.00 

1077.63 

.3323 

.6950 

.0000 

.0003 

• 0300 

.0300 

2. so 

UPS. 99 

.3399 

.5951 

.0000 

."000 

.0000 

.0000 

s.oo 

1393.9C 

.3395 

.6935 

.3333 

.0333 

.3330 

.3300 

10.00 

1S90.9D 

.3160 

.6620 

.0000 

.0033 

.0000 

•sooo 

IS. 33 

175.9.93 

.325S 

.6795 

.3333 

.0033 

.0300 

.0303 

?a.o9 

1673.90 

.3323 

• 668C 

.0003 

.0033 

.0300 

.0000 

2S.33 

1937.53 

.3375 

.6625 

.0000 

.0009 

.0300 

.0000 

3S.UD 

1991.90 

.3955 

.6595 

.0000 

.0033 

.0000 

• oooO 

63.3} 

23CC.0C 

.3503 

.6500 

.0000 

.0000 

.3000 

.0300 

(<S.OO 

1991.90 

.39 55 

.6595 

.3333 

.0333 

.3300 

.3303 

7S.00 

1937.50 

.33>S 

.6625 

.0000 

.0000 

.9000 

.0000 

63.33 

1970.93 

.3323 

.6663 

.3333 

.0331 

.3330 

.3303 

6S.U3 

1759.90 

.3255 

.6795 

.3000 

.0033 

.3000 

.0300 

PJ.33 

1593.93 

.318) 

.6623 

• 3333 

.0339 

.3000 

.3303 

6S.33 

1393.90 

.3395 

.6905 

.ODDO 

.0003 

• 3000 

• DODO 

97.53 

1165.99 

.3399 

.6951 

.0003 

,0003 

.3300 

.0000 

99.00 

1377.63 

.3323 

.6980 

.0003 

.0003 

.3330 

.0000 

99.50 

1039.90 

.3313 

.6990 

.0000 

.0009 

.0390 

.0000 

100.30 

1000.30 

.3333 

,7'’33 

,3333 

.0333 

• 3300 

.3303 


PROFILE 

tepp 

axial 

axial 

WALL 




PiPtPET 

EPS 

E8P0NE91 

TE«P 

PRESS 

TEMP 

RAVENUMSEP 

LAST 

RMS 



AT 

TC.K 

P0.A7M 

79, K 

ll/CMl 

XTERAT109 

RESIDUAL 

P»9A«ETE9S 

PA 

.05000 

1939.37 

.33712 

looo.oo 

3130. 

3 

.996*37 

PA9AWE7ESS 

PA 

.13313 

1928.85 

.39025 

13)3.30 

3133. 

2 

.126*06 

PA9*“ETE9S 

PA 

.15303 

1931.98 

.39395 

1000.00 

3130. 

2 

.899-07 

PABAHETE9S 

PA 

.23000 

1976,28 

.39659 

1009. 00 

3100. 

Z 

.996*37 

PA9A“ET£BS 

PA 

.33003 

2321. IB 

.36001 

10)3.30 

3130. 

3 

.332-)9 

PA9A«ET£0S 

PA 

.93303 

2063.31 

.39PDO 

1003.30 

3130. 

3 

.225-35 

papametebs 

PA 

.53303 

2135.99 

.39832 

1333.30 

3133. 

2 

.767-09 


/ 
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•l*tMEFE9S 

PA 

.63303 

2144.81 

.34*36 

lOOO.OO 

3inj. 

2 

.437-04 



PtRAHCTetS 

PA 

• 0S0C3 

1614. 'i 

• 36718 

1000.00 

3630. 

3 

.722-95 



PltAXeTCRS 

PA 

.13000 

1164,00 

.36065 

1030.00 

3630. 

2 

.118-16 



R*R»**ereRs 

PA 

•tscoo 

I92S.OO 

.36566 

1000.90 

3530. 

2 

.261-96 



RARAMEtERS 

PA 

.23303 

1R6S.TT 

.36223 

1313.00 

3630. 

2 

.144-06 



RARAMEfERS 

PA 

.33003 

2031.56 

.34620 

1000.00 

3600. 

2 

.545-04 



PARAMETERS 

PA 

.49000 

2091. 5S 

.34035 

1000.00 

3690. 

2 

.727-36 



parameters 

• A 

.S3003 

2160.02 

•33688 

1000.30 

3600. 

2 

•617-36 



PARAMETERS 

PA 

.63000 

2205. 4T 

.33533 

1000.30 

3630. 

2 

.257-35 



PARAMETERS 

PA 

.35003 

1639.94 

.3S54S 

1330.30 

3775. 

3 

.271-0‘t 



parameters 

PA 

.13000 

1663.20 

.36384 

1000.00 

3776. 

2 

.230-96 



PARAHETERS 

PA 

.IS 000 

1923.19 

. 36244 

1000.00 

3776. 

2 

.634-97 



parameters 

PA 

.23000 

1962.10 

.35107 

1003.30 

ST76. 

2 

.684-37 



parameters 

PA 

.33003 

2034.16 

.I46S0 

1000.00 

377S. 

2 

.117-94 



parameters 

PA 

.*3903 

2999.36 

.34027 

1133.30 

3776. 

2 

.663>D5 



parameters 

PA 

.53003 

2162.49 

•33S81 

1000.00 

3775. 

2 

.278-96 



parameters 

PA 

.63000 

2222.74 

.33276 

1000.00 

3775. 

2 

.741-96 



PARAMETERS 

PA 

,05000 

1862.26 

.34533 

1033.30 

4190. 

3 

.179-34 



parameters 

PA 

.13009 

1912.02 

.34667 

1000.90 

4130. 

2 

.659-97 



parameters 

PA 

.15009 

1941.17 

.34788 

1313.30 

4130. 

2 

.114-06 



parameters 

PA 

.23303 

1979.48 

.34898 

1333.39 

4130. 

2 

.544-07 



parameters 

PA 

•33903 

2025.86 

•34942 

1000.00 

4100. 

2 

.538-04 



parameters 

PA 

.43000 

2086.02 

.34815 

1000.00 

4100. 

2 

.966-96 



PARAMETERS 

PA 

.59009 

2142.23 

.34760 

1339.30 

4190. 

2 

.648-16 



PARAMETERS 

PA 

.63000 

2196.26 

.34756 

1009.90 

4193. 

2 

.429-06 





OfPEVOENT parameter TC.K WITH 2ETAS - 

.57*03 




NT 

uavenum<:r 


•053 

.too 

.ISO 

.200 

.303 

.400 

.500 .690 

.690 

.600 

IJ/CM) 











3100. 

201 B 

•593 2314, SIT 2009.089 

2304.835 

1992.628 

1977. 6S2 

1962,721 1944.929 

• uOO 

.000 

3539, 

194 9 

•232 1969.663 1982.106 

1992.329 

2333.023 

2935.885 

7997. 253 2005.689 

• 090 

• 000 

3TT6, 

1954 

.161 1865,863 1980.276 

1990.664 

2305.611 

2013.694 

2019,719 2022.863 

.090 

.000 

4100. 

1996 

.476 1997,681 1998.276 

1999.334 

7333.398 

2'333.361 

1999.428 1996.399 

.090 

.000 


dFtN 


page is 
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APPENDIX K 

EXAMPLE OP PROHLr. WITH HtGH-TEMPERATURB CORE 
Data Cards With Meaiurad Data, Two Wave Numbers, ProHle PROF ■ i! 
No Boundary Layer, Gases H2O and N2 


IM 4 WI woier:,p»or:j,iRta 3 :i,paTH:ioc.*iU 2 ,uaVtsS 63 n.,sion., 6 * 3 ». k 

aai:,:3*tl3,,a33a!3,6«3.,HltULStM31b,.hISZ,b*n. , 

t V: l.,Ce«a:3>.,Tu:U 33., TC:|'&('.,NT:. 1.3,1 S3 :U3C.,Jb:S,UPlIH:l., 3330 1, 

1.. 30 nO., 30 UIM:,n 2 .fc 33 ... 0 ?, 60 r.,Utltl?s 2 »P?^il..Pes.! 4 ,.b 6 »a* 3 .. 

N»;7».3iUs3.,ps:, 3« •7,»*),,VK:6iTaTaf>:t33.il333. «tS33.f ?C3S.t?S30. t 

3560.. 2*B.,NSrb,I3Tafc:b03.»ICnC.,16CO.,2000.»?5a0,,»333.*2*O.f 

Ul525,2s3.i.31,.D2,.U5,#l,#2,.3,.a,.b,.6,.?t.3t3al.i 


13l,RblC,T*ir03. « 


H70 

N2 





1. 

.2 





.2^3 

.U36 

.333 

.315 

.238 

' .2b5 


.3157 

.3372 

.Obb'4 

.133 

.127 

.i?b 

.S2b 

1.13 

2.15 

3.35 

a. 36 

.uva 

.113 

.23b 

.bb3 

2. as 

18.0 


Printout oi Input Data 


laiai 

H03L 

- 

♦ 2 





pPSf 

z 

♦ 2 





1P6»3 

z 

♦ 1 





pith 

z 

.13333rc36*33 





II 

z 

♦2 





wlVi 

z 

.35B333C36*3a, 

. 310000036*34, 

.003930306*00, 

.000000006*90, 




.13J333L36433. 

.333333036*33, 

.033133306*30, 

.303039006*90 


HRL 

z 

.31133003- >32, 

.343300006-02. 

.000900306*30, 

.000000006*00, 




.D00a30U3f*30, 

.•300000006*00, 

.000930306*30, 

.300000006*00 


HUUL 

z 

.a31633C3E«33. 

.«352D3036*33, 

.033933336*39. 

.333099996*99, 




.300330r3E«90i 

.000300006*03, 

•000939306*90, 

.000000906*00 


rasa 

z 

. 133333L3C01 





CRPk 

z 

.50003000601 





TW 

z 

.12D93000643a 





TC 

z 

.16b333C.9C*9a 





M 

z 

.600030006*90 





TUB 

z 

.110030036*9a 





JS 

z 

*6 





UPLI»’ 

2 

.100990036*91. 

.300000036*04, 

.103900906*01, 

.300033336*34 


DDTI.IM 

2 

.203930036-31 . 

.600000006*03. 

.203900306-91, 

.500000006*03 


LI 

2 

*l 





L2 

2 

*2 





PS 

2 

.100030036*31 





PC 

2 

.353339L96*33. 

.650303036*3), 

.093339306*30. 

,000000006*00. 




.300930036*90. 

.000000006*00 




NP 

2 

.500030036*90. 

.500000006*00, 

.000900306*00, 

.000000006*00. 




.393993t’96*33. 

.330333036*33 




Ptl 

2 

.303030096*30. 

.700000036*00. 

.000900306*00. 

-000000006*00. 




.333330006*30. 

.000000006*00 




aa 

2 

*b 





Tktap 

2 

.b0003CD36*93( 

.100000036*04, 

.153930306*94, 

.200000306*04. 




.253330C36t3a. 

.300303096*34, 

.033333396*00, 

.000000006*00 


NB 

2 

♦ 6 





TBTAP 

Z 

.bC0030006*33. 

.100000006*04 , 

.150000906*04, 

.200000006*04, 




.253330096*3a. 

.133393036*34, 

.093339)96*30, 

.933093996*09 


J1 

2 

♦ 25 





2 

2 

.000090006*90. 

.100000006-01, 

.203930306-01. 

.500003006-01, 




.100030006*33, 

.200300036*3), 

.333333306*30. 

.433093996*03, 




.503390036*93, 

.600000006*03, 

.703030306*00. 

.900000006*00, 




.193333136*31. 

.130309006*01, 

. 100000306*31 • 

.100000006*01, 




.100030036*31. 

.103030036*31 

.100000006*01, 

.100330336*31. 

.100090096*01, 


IT5P 

2 

♦ 13 





T« 

2 

.300330036*33 





Tcao 







Gtb: 


K20 K2 





4LPHI 


1.03 .20 

.00 ,00 

.30 .03 

.-.0 

.39 

KTABr 


.2563 .aSBC 

.3030 .3150 

.2080 .2653 

.3000 

.0900 



.30*0 .0137 

.0172 .0854 

.1’33 .1273 

.3993 

.0)39 

BT»B: 


.326 .525 

1.100 2. ISO 

3.353 4.353 

.933 

.399 



.374 .113 

,23b .550 

2.453 18.003 

.909 

.900 
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ComputitloD ud ProNIti 


••on It 

TEMP 

SXISL 

SXISL 

USLL 


•ALL PRESS 





•••INEURt 

eXPORERT 

te«p 

PSF&S 

TEMP 

PRESS EXRONCNT 






NT 

»c,« 

PC.STH 

T4.H 

PH.STH 4P 






•soooc 

isoo.no 

.30000 

1 1900.00 


.30000 .00000 





IRCOEMERTS 3Pt 

.OOOCO 

.00 

.00000 

.00 

IIERSII04: 

0 RMS residual: 

.349*30 

STEP: 

.309 


•••meiets •* 

•S30C3 

1109.90 

.30Q0O 

1900,00 

RESIOUSLS: 

-.974*30 -.904.30 

.160-91 

.000 

.000 






nSVENUHOEAl 1/CHI s 

3639. 3103. 

3500. 


D. 

P. 

INCDEMERTS }P« 

. 10903 

«?.»6 

•.OOOOT 

190.00 

nCRSTtONs 

1 RMS residual: 

.920*90 

SICP: 

.193*00 


•••tMEFEftS Pi 

.«3000 

11*3. 06 

.34643 

1320.00 

residuals: 

-.176*33 -.139*33 

• 192-9I 

.009 

•000 






4SVENUHBERII/CHIS 

3000. 3109. 

3000. 


0. 

0. 

INCRFMENTS 9P« 

.00003 

06.10 

•.C0019 

139.00 

ITERSTI04: 

9 RH$ residual: 

.646-01 

S7FPS 

.133*00 


POttMETCtS P» 

.00009 

I9«B.B| 

.34990 

1409.00 

RESIDUALS: 

-.S13-0I -.399-91 

.231-02 

.ODD 

.000 






nAVeNUHBERtl/CHi: 

3S39. 3109. 

3500. 


0. 

r. 

INCRFMEtiTS OPS 

.93003 

94.63 

.00066 

47.19 

ITERSIIor: 

7 RHS residual: 

.969-03 

STEPS 

.455-01 


P«9*H£TePS PS 

.03009 

1673. TO 

.34991 

1499,19 

RESIOUSLSs 

.747.93 .636.33 

-.198-03 

.000 

.000 






4SVC HUHBeRtl/CMi; 

3000. 3100. 

3500. 


0. 

0. 

INCREHERTS ops 

.33303 

• •49 

.00010 

-.09 

itersiion: 

4 RMj residual: 

.999-07 

STEPS 

.060-0 3 


P*»S«EFE#S PS 

.03009 

1973.99 

.30000 

1490.67 

RESIDUALS: 

.596-07 .447-97 ' 

-.605-07 

.000 

.000 






XSVEHUHBERII/CMIs 

3039. 3103. 

3500. 


0. 

0. 

PSRFRCrCP E9ROR 




USVERJ4BE4I1/CHI PESIDUSL, PERCENT 





PEftCENF 

.91 

•.71 

.91 

3. St 

3530. 

1.6t 






.3 

9.6 

.0 

-T.O 

3133. 

3.3 






.3 

•.1 

6.3 

. 6 

3SOO. 

1.4 






.3 

.3 

.0 

.0 

0. 

.0 






.0 

.0 

.0 

.0 

Q. 

.0 





QMS rorsL 

.3 

9.9 

6.3 

7.6 

0. 

.0 





PSVENUHBER(1/CH| 

= 3033. 

3133. 

3. 

3. 

3. 3 

. 3. 0. 





tpsnsmittsnce 

: .«32 

.630 

.000 

.030 

.030 .000 .000 .000 






PROFILE 


PFRim PRCS&ORE PROFILES E»tH SFSftTH 


X.CM 

fEMP.K 

490 

49 





.00 

1100.00 

.3300 

.7000 

.0000 

.0000 

• 0000 

• 0000 

.50 

1166.19 

.3313 

.6990 

.0333 

.093) 

.3390 

.9339 

1.00 

1269.39 

.3399 

,6900 

.0333 

.0933 

•0300 

.0303 

9.00 

1445.96 

.3349 

.6901 

.0003 

.0039 

.0000 

.Q900 

0.33 

1066.60 

.3390 

.6905 

.0333 

.0039 

.0000 

.0000 

10.00 

1669.03 

.3180 

.6890 

.0000 

.0000 

.0000 

.0000 

10.30 

1740.73 

.3905 

.6740 

.0000 

• 0000 

.0390 

.0000 

20.00 

1809.4 3 

.3399 

.6660 

.3333 

.0999 

.0900 

• 3303 

25.00 

1604.63 

,3575 

.6620 

• 0000 

.0009 

.0000 

.0000 

30.00 

1897.35 

.3493 

.6060 

.0333 

.0939 

.3300 

.9)03 

30.00 

1930.57 

• 345S 

.6545 

.0000 

.0099 

• 3000 

.0000 

43 .9] 

1904.39 

.3483 

.6520 

.0000 

.0000 

.0300 

.0000 

50.00 

1973.99 

.3SC0 

.6000 

.0000 

.0009 

.0300 

.0000 

69.33 

1904.39 

.3463 

.6020 

.0003 

.0000 

.0300 

.0000 

66.00 

195C.57 

.3450 

.6540 

.0000 

.0033 

.0300 

.3309 

70.00 

1697.30 

.3499 

.6560 

.0000 

.0000 

• 0000 

.0000 

76.00 

1804.63 

.3375 

.6625 

.9333 

.0333 

.0300 

• 3909 

60.00 

1639.43 

.1390 

.6660 

.0000 

.0090 

.0000 

.0000 

66.00 

1740,73 

.3900 

.6745 

.0333 

.0033 

.0300 

>0909 

90.00 

1669.53 

.3169 

.6820 

.0003 

.0033 

.0000 

.0000 

96.33 

1566.85 

.3396 

• 6930 

.ODOD 

.0000 

•oono 

.0000 

97,50 

1440,96 

,3049 

.6951 

.0000 

.0009 

,0000 

.0000 

99.33 

1969.39 

.3393 

.6960 

.0000 

.0000 

.0300 

.0000 

99.00 

1166.19 

.3313 

.6990 

.0000 

.0939 

.9300 

.0309 

log. 00 

1100.00 

.3303 

.7000 

.0000 

.0000 

• 0300 

• 0000 


OIOfilNiU* PAGE! IS 
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i ABLE I. - COEFFICIENTS TO BE USED IN EQUATION (16) 



Al 

^2 


A 4 

^5 

0 

- 0.68 

-0. 35 

0. 55 

0.07 

1.00 

2 

-.37 

.62 

.71 

.43 

1.37 

5 

.60 

3.15 

.90 

.79 

1.93 

10 

3.19 

9.25 

1. 15 

1.30 

2.85 

20 

9.10 

22.8 

1.63 

2.26 

4. 70 

50 

29.1 

67.1 

2.78 

4.56 

10. 25 

100 

105 

209 

4.05 

7.09 

19. 50 


TABLE n. > SPECTRAL UNE BROADENING COEFFICIENT, 

a* (FROM REF. 16) 


Absorbing gas 

Broadening gas 

HgO 

COg 

CO 


N2 

tSM 

°2 

HjO 

1 

0.77 



0.20 

QQI 

1 

COg 


1 



.77 

1.17 

0.81 

CO 



1 

1.12 

.95 

.85 


CH 4 


1.25 


1 

.75 
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TABLE III. - ESTIMATED PERCENT ERROR OF T„ AND T* 

c w 

AND TRANSMITTANCE FOR VARIOUS PATH LENGTHS 
FOR PROFILE OF FIGURE 1(b) 

« 

/ [Measurement wavenumbers 3100, 3400, 3500, and 4100 cm* 

Tg = 2000 K, n,j. = 0. 5, = 0. 35 atm, = 0. 30 atm, 

Up s 0. 5. Errors caused by one percent error of radiance 
measurement, and five percent error of tabulated absorp- 
tion coefficients^ 


j Percent errrr in T^, 

msH 

Path length, cm | 

K 






50 

100 

200 

400 

1200 

3.0 

2.1 

1.4 

1.0 

2000 



11.5 

^ 

6.3 

4.2 

3.2 

Percent error in T^ 

1200 

14.1 

9.7 

n 

6.5 

4.4 

2000 

20.0 

9.8 

5.7 

3.8 

Transmittance at 3500 cm"^ 

1200 

0.703 

0.591 

0. 465 

0.333 

2000 

.713 

.602 

.476 i 

.343 
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TABLE IV. - PARAMETER INDEX CONTROL 


KEY 

KEY1(KEY) 

KEY2(KEY) 

KEY3(KEY) 

Variable parameters 


kx 

K2 

K3 

Fig. 1(a), 

Fig. 1(b), 





PROF:=l 

PROP=2 

1 

1 

3 

1 

“t 

Te Pc K 

2 

1 

2 

1 

Ht Tg 

Pc 

3 

1 

3 

2 

Ht Pq 

T T* 

C 

4 

1 

1 

1 


Tc 

3 

2 

3 

1 

TcPc 

Pct; 

6 

2 

2 

1 

Tc 

Pc 

7 

3 

3 

1 

Pc 

T* 

W 

8 

0 

0 

0 

None 

None 


NASA'Law it -Cotn'l 



0 


NOn-lsothermal temperature profiles, 
Tyi, ■ constant. 

.5f— 



(&) Near-Isothermal temperature profiles, with 
Tyj, • constant and boundary layer thickness 
Xj • constant. 



(cl Partial pressure profile for each absorbing 
gas. 


Figure 1. - Assumed profile functions. 


/ ‘-3500, 3550 
1775 


.2 


.3 

Exponent, nj 


.5 


Figure 3. - Test of wavenumbers for profile of figure 1(a) with 1 
1000 K. Pp ■ 0. 35 atm, Py^ • 0. 30 atm. n,, - 0. 50, path length 
lOO cm Hydrogen-air combustion at 1 Itm. 





Figure 4. ■ Test of wavenumbers for profile of figure 10) with • 
0. 5, Pc ' 0, 35 atm, Pj^ > 0, 30 atm, np • 0. 50. path length 100 cm. 
Hydrogen-air combustion at 1 atm. 



C 


ia 



Figure 5 . - Estimated error in due solely to errors in radiometer 
and tabulated absorption coefficients. Profile of figure 1(a), with 
Tj • aoo K. Tm • 1000 K. Pg • 0. 35 atm, Py. • 0. 65 atm, ‘ 0, 50. 
path length 100 cm, Hydrogen-atr combustion at 1 atm. Error of 
radlometry 1 percent, error of tabulated absorption coefficient 
5 percent. 
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I 


ETForin T, 



Figure 6. *■ Estimated error In T*. due solely to errors In radlometry 
and In tabulated absorption ojefflclents. Profile of figure 1(b) with 
Tc • 2000 K, nr - 0. 5, p. • 0. 35 atm. p^ • 0. 30 atm. iTn - 0.m _ 
path length llw cm. Hydrogen-air comoustlon at 1 atm. Error of 
radlometry 1 percent, tabu lated absorption coefficient 3 percent. 
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incremented 

paremeters 


KEY«1 


Figure 7. - Flow chart for PROFIL. 
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Figure 8. » Flow chrrt for WAVES. 










» 







(c) Profile of figure 1(b) with boundary layer, 
n-f = 0. 5. = 1100 K. 

Figure 9. - Example of profile with high- 
temperature core. Wavenumbers of 

measurement 3500 and 3100 cm . ORIGINAL PAGE IS 

OP POOR QUALITY 




